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Abstract: 

We develop a general theoretical framework for the description of higher-twist baryon 

operators which makes maximal use of the conformal symmetry of the QCD Lagrangian. 

The conformal operator basis is constructed for all twists. The complete analysis of 

the one-loop renormalization of twist-4 operators is given. The evolution equation for 

three-quark operators of the same chirality turns out to be completely integrable. The 

^ spectrum of anomalous dimensions coincides in this case with the energy spectrum of the 

twist-4 subsector of the SU(2, 2) Heisenberg spin chain. The results are applied to give 

a general classification and calculate the scale dependence of subleading twist-4 nucleon 

distribution amplitudes that are relevant for hard exclusive reactions involving a helicity 

flip. In particular we find an all-order expression (in conformal spin) for the contributions 

O of geometric twist-3 operators to the (light-cone) twist-4 nucleon distribution amplitudes, 

OO 
O 

> 

•i-H 



which are usually referred to as Wandzura-Wilczek terms. 
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1 Introduction 



Higher-twist effects in hard processes in QCD generically correspond to corrections to 
physical observables that are suppressed by powers of the hard scale. They are important 
in order to achieve high accuracy, and interesting because higher-twist corrections are 
sensitive to fine details of the hadron structure. A theoretical description of higher 
twist effects within QCD factorization involves contributions of a large number of local 
operators which are much more numerous compared to the leading twist so that the 
choice of a proper operator basis is important. This choice is not unique, as exemplified 
by the two existing classical approaches to the twist-4 effects in deep-inelastic lepton 
hadron scattering [TJ [2] . The "transverse" basis of Ref . j2] leads to simpler coefficient 
functions whereas the "longitudinal" basis of Ref. [TJ (see also [3]) allows for a parton- 
model-like interpretation [I]. 

The renormalization of higher-twist operators corresponds to the scale dependence 
of the physical observables. For twist three, the corresponding study is essentially com- 
pleted. The anomalous dimension matrix for baryon operators was first calculated in [5], 
for chiral-even quark-antiquark-gluon and three-gluon operators in [6], and for chiral- 
odd in [7]. The structure of the spectrum of twist-three anomalous dimensions is well 
understood |8j [9j QUI El EEl [El El 02] and in some cases explicit WKB-type expansions 
are available that allow to calculate anomalous dimensions to arbitrary accuracy: The 
size of the mixing matrix plays the role of the expansion parameter. Beyond twist three 
much less is known. Up to now, anomalous dimensions have only been calculated for 
a few operators of lowest dimension (e.g. [151 E32])- m addition, the structure of the 
most singular parts of the mixing kernels for small values of the Bjorken variable that 
are relevant for the contribution of two-pomeron cuts in high-energy scattering processes 
was considered in [HI [19] . 

The modern approach for the calculation of leading-order anomalous dimensions of 
higher-twist operators makes maximal use of the conformal symmetry of the renormaliza- 
tion group equations. Historically, the importance of conformal symmetry in the present 
context was first understood for the leading twist pion distribution amplitude and it was 
instrumental for the proof of QCD factorization for the pion form factor [201 |2Tl [22]. 
A general formalism was developed in [23] for the special class of so-called quasipar- 
tonic operators that are built of "plus" components of quark and gluon fields. For each 
twist, the set of quasipartonic operators is closed under renormalization and the renor- 
malization group (RG) equation can be written in a Hamiltonian form that involves 
two-particle kernels given in terms of two-particle Casimir operators of the collinear 
subgroup 5X(2,R) of the conformal group. In this formulation symmetries of the RG 
equations become explicit. Moreover, for a few important cases the corresponding three- 
particle quantum-mechanical problem turns out to be completely integrable and in fact 
equivalent to a specific Heisenberg spin chain [8j. An almost complete understanding 
achieved at present of the renormalization of twist-three operators is due to all these 
formal developments, see [2"H |2"5] for a review and further references. 



2 



The goal of this paper is to generalize some of the above techniques to the situation 
where not all contributing operators are quasipartonic, as it proves to be the case starting 
with twist four. Apart of the needs of practical applications to QCD phenomenology, 
our work is fuelled by the recent study [261 [27] where it was shown that diagonal part of 
one-loop QCD RG equations (for arbitrary twist) can be written in a Hamiltonian form 
in terms of quadratic Casimir operators of the full conformal group 50(4, 2) instead of 
its collinear subgroup. Moreover, all kernels can be obtained from the known kernels for 
the collinear SL(2,M) subgroup [23J. Although much of the formalism appears to be 
general, in this paper we concentrate on the simplest example of non-quasipartonic twist- 
four baryon operators that contain two "plus" and one "minus" quark field, schematically 

q+q-q+ , 

and their mixing with (quasipartonic) four-particle operators involving a gluon field, of 
the type 

Our main results can be summarized as follows. 

First, we construct a complete conformal operator basis for arbitrary twist, with 
"good" transformation properties. We then specialize to the case of twist-4 baryonic op- 
erators, calculate all one-loop evolution kernels including the mixing with four-particle 
operators involving a gluon field, and check that the kernels are SL(2) invariant, as ex- 
pected. The operators involving three quark fields with the same chirality do not mix 
with the operators involving both chiral and antichiral quarks, so that these two cases can 
be considered separately. The evolution equation for three-quark operators of the same 
chirality turns out to be completely integrable. The spectrum of anomalous dimensions 
coincides in this case with the energy spectrum of the twist-4 subsector of the SU(2, 2) 
Heisenberg spin chain, confirming the prediction of [27J. For both cases, we present a 
detailed study of the spectra of the anomalous dimensions. Finally, these results are 
applied to give a general classification and calculate the scale dependence of subleading 
twist-4 nucleon distribution amplitudes that are relevant for hard exclusive reactions 
involving a helicity flip. In particular we introduce novel four-particle distribution am- 
plitudes involving a gluon field, and derive explicit expressions for the expansion of all 
distribution amplitudes in contributions of multiplicatively renormalizable operators in 
first three orders of the conformal expansion. As a byproduct of our analysis, we give an 
expression for the contributions of geometric twist-3 operators to the (light-cone) twist- 
4 nucleon distribution amplitudes, which are usually referred to as Wandzura-Wilczek 
terms. 

The presentation is organized as follows. We begin in Sect. 2 with a short exposition 
of the spinor formalism that is used throughout our work. This formalism is standard 
in the studies of SUSY theories but is used rarely by the QCD community so we felt 
that a short summary is necessary. Next, conformal transformation properties of the 
fields are considered in some detail. A complete basis of one-particle light-ray operators 
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is constructed for chiral quark and self-dual gluon fields in QCD, cf. Eq. (2.59), which 
is one of our main results. In Sect. 3 we specialize to the particular case of baryonic 
twist-4 operators which are the main subject of the rest of the paper. Renormalization 
group equations for the light-ray baryonic operators are derived in Sect. 4. We discuss 
general properties of these equations, give a summary of the relevant conformal invariant 
evolution kernels, introduce a convenient scalar product on the space of the solutions and, 
finally, give explicit expressions of the Hamiltonians for all cases of interest. Solutions 
of the renormalization group equations for twist-4 operators are considered in Sect. 5. 
For three-quark operators of the same chirality the problem turns out to be completely 
integrable. We find the corresponding conserved charge and discuss the relation of this 
result to the approach of [26| [27] . A simple analytic expression is found for the lowest 
anomalous dimension in the spectrum of chiral quark twist-4 operators with odd number 
N = 2k + 1 of covariant derivatives. For other cases the spectra are studied numerically. 
The results are presented in the Figures and for the first few iV also in table form. 
It turns out that differences between twist-4 and twist-3 operators mostly affect a few 
lowest eigenstates (for a given N); the upper part of the spectrum is universal: the 
anomalous dimensions appear to be almost independent on twist and chirality. Explicit 
expressions for the nucleon distribution amplitudes taking into account first three orders 
in conformal spin and Wandzura-Wilczek corrections are given in Sect. 6. The final 
Sect. 7 is reserved for summary and conclusions. 



2 Spinors and Conformal Symmetry 

For applications it is important to have an operator basis with good transformation 
properties with respect to the collinear SL(2,M) subgroup of the conformal group. It 
is well known that analysis of tensor properties of operators is greatly simplified in 
the spinor representation. Although this formalism is standard, a number of different 
prescriptions exist in the literature for raising and lowering indices, normalization etc. 
In order to make our presentation self-contained we choose to begin with a summary of 
the definitions and basic relations of the spinor algebra, and also introduce some general 
notation that is used throughout the paper. Our conventions are similar but not identical 
to the ones accepted in Ref. [28] . 



2.1 Spinor formalism 

The Lorentz group 5*0(3, 1) is locally isomorphic to the group of complex unimodular 
2x2 matrices, SL(2, C). To make this explicit, each covariant four- vector x can ue 
mapped to a hermitian matrix x 

( x + x 3 xi - ix 2 \ _ u (n . x 

x = = ay 7 ( 21 ) 
\x 1 + %x 2 x - x 3 J 
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where cr M = (1, a) and a are the usual Pauli matrices. A Lorentz transformation 
x' = A fJi u x u corresponds to a rotation x' = MxM\ where M 6 SX(2,C), and the 
homomorphism A — > M defines a two-dimensional (spinor) representation of Lorentz 
group, u' = Mu. The correspondence between A and M is not unique and in general one 
might consider four representations defined by the homomorphisms A — > M, M*, M~ 1,T 
and M~ x ^ . The vectors from the corresponding representation spaces - spinors - are usu- 
ally denoted and u a , respectively, i.e. u' a = Mjup, v! h = Mfup etc. The 
representations M and M~ 1,T (also M* and M~ 1 ^) are equivalent since a 2 M = M~ 1,T a 2 . 
The intertwining operator a 2 is proportional to the Levi-Civita tensor e. We define 



.12 



— e 



.12 



ei2 — e — ~ e i2 

and accept the following rule for raising and lowering of spinor indices (cf. [28J) 



(2.2) 



u 



u, 



u p tp a , 



Ur 



which is consistent with (2.2). Note that e a ^ = —e^ a = 5f and e a k = —£n a = 5°t 



(2.3) 



When it is not displayed explicitly it is implied that undotted indices are contracted 



'up-down' : 



{UV} 



dcf 



-u a v a and dotted ones "down-up" , (uv) = f UaV a 



Next, we define (u a )* = and («")* = u a that is, again, consistent with (2.2) and 
results in (uv)* = (vu). The Fierz transformation for Weyl spinors reads 



(u x u 2 )(viv 2 ) = (uiv 1 )(u 2 v 2 ) - (uiv 2 )(u 2 v 1 ) 
which is a consequence of the identity 

^ab^cd ^ac^bd ^ad^bc 



(2.4) 



(2.5) 



In addition to cr^ = (1, a) it is convenient to introduce (a M ) a/3 = (1, — a) so that 



(p^a-P — (a' M ) l3a , and define x = x^, cf. (2.1). One easily finds that 



-a b° 



For completeness, we give below some useful identities involving matrices: 



2g 



(2.6) 



( a »<j» + cT u a%P=2g^5 a , 
Generators of the Lorentz group read 



( a ) p= ° ~ J /3' 



(2.7) 



(2i 



5 



or, in the explicit form 

„0i • i Jik -Oi • i -ik Aj-ki„i (n n\ 

a = —ia , a = le J a J , a = in , a = le J cr J , 

They satisfy the self-duality relations 

a ,u = l _ e ^ 0pw ; ^ = - % -e^a^ . (2.10) 

where e i23 = 1- 

A four-dimensional Dirac bispinor is written as 

and the 7 M matrices take the form 



q = i^M (2.11) 



For the common a^ u = |[7 M ,7^], 75 = i7°7 1 7 2 7 3 and the charge conjugation matrix 
C = ry 2 7° one finds 



a 



mv- -* = {° *})• c= \ -<*».» ■ (2 - i3) 

Irreducible representations of the Lorentz group are labeled by two spins (s, s). The 
representation space is spanned by tensors T ai a2 h i which are symmetric in dotted 
and undotted indices separately. In particularly, the Weyl spinors ip (chiral) and \ 
(antichiral) belong to the representations (1/2,0) and (0, 1/2), respectively, whereas the 
Dirac spinor transforms as (1/2, 0) © (0, 1/2). 

The gluon strength tensor transforms as (1,0) © (0, 1) and can be decomposed 

as 

F afS,a/3 = a aa a p,3 F ^ = 2 { € apfafi ~ Ca/j/d/j) ( 2 - 14 ) 

where f a p and f & p are chiral and antichiral symmetric tensors, /* = /, which belong to 
the representations (1,0) and (0, 1), respectively. One obtains 

u=l<^„, /<* = 4*2^ ' ■ (2 - 15) 

or in terms of the gauge field A aa : 

U = \ {D a * A*, + DfA aa ) , f. $ = \ (D a a A a$ + D«A aa ) , (2.16) 
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where the covariant derivative is defined as = <9 M — igA^. The expressions for F^ v 
and the dual strength tensor F^ u = -^e^ upa F p(T are 

= \ (<T P - <';r ') , ^ = \ (<;r* + • (2.17) 

The Dirac equation for the quark fields reads 

D^ a {x)=Q, D a&x & (x) = (2.18) 

where the covariant derivative is defined as D p = — igA^. The equation of motion 
(EOM) for the fields /, / becomes 

Dffo = g feT<% + XP T a Xp) , Dffy = g feT^ + Xp T a X$ ) . (2.19) 

The class of the operators which are proportional to the equation of motion is closed 
under renormalization (for a more precise statement see e.g. Ref. |29j). On-shell matrix 
elements of such operators vanish and one can consider two operators which difference 
is an EOM operator as being equivalent. 
The equation 

T ai ...a n A-A = a Z$i ' ' ' al Z$ n T ^-^ ( 2 ' 2 °) 

establishes the relation between generic tensors in the usual vector and spinor representa- 
tions. The symmetrization over spinor indices is most conveniently achieved contracting 
the open indices with an auxiliary spinor £. We define 

r c = . . . T ai ___ anA __ A f 1 . . . P . (2.21) 

In particular 

x 5 = m = & = h tec , (2.22) 

etc. 

It is obvious that a symmetric tensor T ai an * * can unambiguously be restored 
from the convolution T% by applying multiple derivatives over £. We define 

= = ep a = 5% , 5% = -t-h = e\ = £ , (2.23) 

so that 

■l) fi d d d d 

Note that the rule for raising and lowering of indices for derivatives over spinor variables 



T ai ...a n ,0i...$ n — „| = | atry, a^n. o^fl. ••' *7h_ ^ ' (2.24) 



is different from that for the spinors themselves, cf. Eq. (2.3): 



— -c — 9 - 9 (9 9^ 
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2.2 Conformal symmetry 

It is known that the QCD enjoys conformal symmetry at the classical level. Although 
this symmetry is broken in the full quantum theory, it leads to strong constraints on 
the form of (one-loop) operator counterterms and will be quite useful in the subsequent 
analysis. The action on the generators of the conformal group on the fundamental fields 
in the spinor representation, $ = ($^,$^) with $g = {V'fjXfj/f} and = /c}, 
takes the form 



z[D,$(» 
i[M a/3 ,$(x) 

z[K ad ,$(x) 



d 



2^ 



$(x 
d 



- d - d 



<9 77 + 2tx adl + 2£ a x 







+ 2£*x c 



(2.26a) 

iD$(x), (2.26b) 

iM^(x), (2.26c) 

<$>(x) = iM^(x) , (2.26d) 
d ' 



^ 1 ~ SQ ~ Q/3 «9^ 



iK a ^{x) , 
(2.26e) 



where <9 ac i = c^d^ and i = 1 is the geometric twist [3T]: for the field with canonical 
scaling dimension £ can and Lorentz spin (s,s) it is defined as t = £ can — s — s. Note 
that we use boldface letters for the generators acting on quantum fields to distinguish 
them from the corresponding differential operators acting on the field coordinates. The 
transformations of the gauge field = A a a£, a !; a are given by the same expressions with 
t = 0. 

In the applications of QCD to high-energy scattering the separation of transverse and 
longitudinal degrees of freedom proves to be essential. It is conveniently achieved by the 
introduction of two independent light-like vectors 



which we choose to be normalized to 

(/iA) = -(A/i) = l, 

Without loss of generality one can take 

A a = (l,0), 
/i a = (0,l), 



n A = 
n 2 = 



n • n 



1/2 



(2.27) 



(2.28) 



= (0,1), 

:-i,o). 



(2.29) 
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Then, for example 



d u = 2(n ■ d) , d 11 = 2(n • d) (2.30) 

are the derivatives in the two chosen light-like directions whereas the remaining two, d 12 
and <9 21 , are the derivatives in the transverse plane. 

Fast moving hadrons can be viewed as a collection of partons that move in the same 
direction, say n M . Whenever this picture applies, quantum fields "living" on the light 

ray 

-> $(zn) (2.31) 

play a special role. Such light-ray fields can be viewed as generating functions for local 
operators that arise through the (formal) Taylor expansion 

= Hzn) = y ^d) k m = ^(d 22 ) k m . (2.32) 

k k 



Note that all local operators on the r.h.s. of (2.32) have the same collinear twist as the 
field $ itself since each d 22 derivative adds one unit of dimension and spin projection, 
simultaneously. We will use a shorthand notation for $>(nz) in what follows. 

With the restriction to light-ray operators the four-dimensional conformal trans- 
formations are reduced to the collinear subgroup SL(2, R) corresponding to projective 
(Mobius) transformations of the line x = zn: 

ab — cd = 1 , 



cz + d 



where a, 6, c, d are real numbers. The generators of the collinear subg roup, S±, Sq can be 
chosen as 

S+=faKji), S- = - % - ( A P A) , So = % - (D - ^ X^M aP - M^tf) , (2.33) 
or, using the convention in Eq. (2.29), 

S+= l -K 22 = l -K'\ S„ = - l -P 22 = - l -P li: S = % -{d - M 21 - M i2 ) . (2.34) 
The explicit expressions are 

S+ =-x 2 ^x l2 d^ + x 22 (t + £ p -^ + P —^J - x^f— - ^2/3?^ 7 (2.35a) 
S- = - l -d 22 , (2.35b) 
So =\ (x 22 d 22 + \ (x 2i d 2 ^ + x 12 d' 2 ) + t + + ?-V) ■ (2.35c) 
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They obey the standard commutation relations 



[S+,S-] = 2S , [S ,S ± ] = ±S ± . (2.36) 

In addition, there exist two operators that commute with all SX(2,R) generators: 

E=i[D+ M 21 + M a ) = i^S" + - (x 2i d 2i + Xli d" + 2t) + f 2 ^ + f^j • < 2 - 37 ) 

=i(M 15 - m 21 ) = i (x 2l a- - x l5 a- + - ^ A _ ^ + pjL j . (2 . 38) 

-E is usually called the collinear twist operator: collinear twist counts the dimension of 
the field minus spin projection, as opposed to the geometric twist t which is dimension 
minus spin. In a slight abuse of language we will refer to H as the helicity operator; the 
name can be justified by observing that for "good" components of the fields (see below) 
the eigenvalue of H coincides with helicity of the corresponding one-particle state. 

A light-ray operator with definite collinear twist E transforms according to the irre- 
ducible representation of the SX(2,R) group with the conformal spin 

j = £ can -E/2. (2.39) 

In particular the SL(2) generators acquire their canonical form 

S + = z 2 d z + 2jz, S = zd z + 3} S- = -d z , (2.40) 

i.e. first order differential operators acting on functions of the light-cone coordinate z. 
The finite form of the group transformations is 

For example, a chiral field ip should be decomposed as 

if>(z) = \i/>.(z) - (Jiil> + (z) , (2.42) 

where 

Mz) = \ a 1p a (z) EE ^(Z) , m + ](z) = ^+(z) , m+]( Z ) = ^+{z) , 

= » a Mz) = M*) , W-](z) = 2^(z) , m_]{z) = - l -^-{z) . (2.43) 



Note that iMi\ appearing in (2.37), (2.38) counts the difference in the number of "first" 
and "second" spinor indices, which is nothing but the Lorentz spin projection on the 
light-ray direction. In particular ip + and correspond to spin projections +1/2 and 
—1/2, respectively. Using explicit expressions in Eq. (2.26) it is easy to check that the 
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^+ 




X+ 


X- 


/++ 


/+- 


/- 


3 


1 


1/2 


1 


1/2 


3/2 


1 


1/2 


E 


1 


2 


1 


2 


1 


2 


3 


H 


1/2 


-1/2 


-1/2 


1/2 


1 





-1 



Table 1: The SL(2,M) spin and twist for the fundamental fields 



fields ip + and ip- indeed transform according to Eq. (2.40) with the conformal spin j = 1 
and j = 1/2, respectively. 

Similarly, for the anti-chiral field x we define the "plus" and "minus" projections as 

X+ = XaX a , X- = X*t (2.44) 

and for the self-dual vector field f a p 

f ++ (z) =\ a \?f a p(z) , /+_(*) =AV , /— (*) =^fap(z) . (2.45) 

The projections for the conjugate fields are defined as ?/>± = (V>±)* etc. The SX(2,IR) 
quantum numbers of the fundamental fields — conformal spin, (collinear) twist and 
helicity — are collected in Table [Tj 

The plus components of the fields, &+(z) = Q\(zn) = {ip + (z),x+{z), f++(z)}, and 
their anti-chiral counterparts — "good" components in conventional terminology — have 
the lowest twist. The product of the plus fields taken at the different points on the light- 
ray 

$ + (2i)$+(2;2)--- < l > +(^iv) 

serves as a generating function^for the so-called quasipartonic operators [B]. An operator 
constructed from N "plus" fields has collinear twist E equal to iV which is the lowest 
possible twist for iV— particle operators. The set of N— particle quasipartonic operators is 
closed under renormalization at the one-loop level. The renormalization group equation 
can be reinterpreted as a Schrodinger equation where the scale \x plays the role of time. 
The corresponding Hamiltonian contains pairwise interactions only and can be written 
in terms of the two-particle Casimir operators of the collinear conformal group [6j. 

The light-ray iV— particle operators containing minus components of the fields, 
{ip-,X-, /+-, . . .}, have twist larger than N and provide one with examples of operators 
that are not quasipartonic. Renormalization of non-quasipartonic operators in QCD has 
never been studied systematically, to the best of our knowledge. On this way, there are 
two complications. 

First, the number of fields ("particles") is not conserved. To one-loop accuracy, the 
mixing matrix of operators with a given twist E has a block-triangular structure as the 
operators with less fields can mix with ones containing more fields but not vice versa. 
Operators with the maximum possible number of fields N = E are quasipartonic. 



"For the moment we ignore the color structure and all issues related to gauge invariance. 
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Second, operators involving minus field components can mix with operators of the 
same twist containing minus, d 11 , or transverse, d l2 ,d 21 , derivatives. These operators, 
therefore, also must be included. The problem is that transverse derivatives generally 
do not have good transformation properties with respect to the SX(2,R) group. In 
concrete applications it may be possible to get rid of such operators using EOM and 
exploiting the specific structure of the matrix elements of interest, e.g. if there is no 
transverse momentum transfer between the initial and the final state. Two well known 
examples are the twist-four contributions to the deep-inelastic scattering (DIS) [TJ and 
to meson distributions amplitudes [32j E3]. The main problem as far as the operator 
renormalization is concerned is that after this reduction conformal symmetry becomes 
obscured. This procedure is also not universal and probably cannot be applied beyond 
twist four. 

In this work we suggest a different, general approach based on the construction of a 
complete conformal operator basis for all twists. In this basis, the SL(2, M) symmetry 
of the renormalization group equations is manifest. To begin with, we will explain our 
construction on the example of a free chiral field ip, the extension to the other fields is 
straightforward. 

Let us examine the action of the SL(2,M) generators in Eq. (2.40) on the light-ray 
operator with a transverse or "minus" derivative, [d 12 ip±](z), [d 21 ip±] (z) and [<9 1: V±](<2). 
It is easy to see that So and S- retain their form, and complications only arise in the 
case of S + which is related to special conformal transformations: 



i[K a6l 1pp]{x) = (x^X^a <9 77 + Ax a6 ) tpp{x) - 2Xf3a1pa(x) 



(2.46) 



cf (2.26). In particular 

i[K 22 ip-)(x) = (x 2 jX j2 <9 77 + 2x 2 2) ip-(x) , 



(2.47) 



The action of the "spin- up" generator S + = iK 22 /2 on the light-ray operator with a 



transverse derivative follows readily from Eq. (2.47) observing that, e.g. 



K 22 ip±\ 



x , 



K 22 , [<9 2 V ± ]] (z) = (> 

and taking into account that & )ia x^ = 25^5^ and x 22 = z. One obtains 

S + [d 2i i; + ] (z) = (z 2 d z + 3z)[d 2i i; + ] (z) , 
S+P^+Kz) = (z 2 d z + 2z)[d li ^ + ](z) , 
S + [d 12 ^ + }(z) = (z 2 d z + 3z)[d 12 ^ + }(z) - 2^(z) 



(2.48) 
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and 



S + [d 2i ^](z) = (z 2 d z + 2z)[d 2i ^} (z) 
S + [d li ^](z) = (z 2 d z + z)[d li ^](z), 
S + [d 12 ^}(z) = (z 2 d z + 2z)[d 12 tfj_}(z) 



(2.49) 



We see that the generator S + take the standard from (2.40) for all cases exce pt fo r 
[d 12 ip + ](z). Fortunately, this derivative can be eliminated with the help of EOM (2.18): 



d 22 ^ 



(2.50) 



The first equation in (2.50) allows to replace all occurrences of d tjj + (z) by —dr i])_{z) 



—d z i[)-(z). The second one, in principle, can be used in either direction since d 21 ip_ and 
d 11 ijj + both have "good" transformation properties. It turns out, however, that eliminat- 
ing d 21 tp- in favor of d u ip + is advantageous since it leads to a simpler complete operator 
basis in a general situation and we adopt this option for what follows. The remaining 
four independent operators d u ip + , d llr ip_, d 2l ip + and d 12 ifi_ transform according to the 
irreducible representations of the collinear conformal group with spin j = 1, j = 1/2, 
j = 3/2 and j = 1, respectively. Note that the "minus" derivative does not change the 
conformal spin of the light-ray operator, whereas a "good" transverse derivative increases 
the spin by 1/2. 

The above construction can be generalized for an arbitrary number of derivatives. It 
is easy to verify that the following fields 



^ m \ z ) =[{d 2i ) 2 i- 2 {d 1{ ) 2m 4> + }{z) 



(2.51) 



transform according to the spin-j representation of the SX(2,R) group, Eq. (2.40). All 



other combinations of derivatives can be reduced to this basis with the help of E OM. In 
particular, all pairs d 12 d 21 can be replaced by d 22 d u which is a consequence of (2.50). 



Next, we consider which modifications have to be done in the (interacting) gauge 
theory. In the first place we have to replace ordinary derivatives by the covariant ones. 



This is achieved by modifying the definition of the light-ray operator (2.32) to include 
the factor 



where 



[0, z] = Pexp 



igz / duA (uz) 



(2.52) 



(2.53) 
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is the light-like Wilson line in the appropriate (fundamental or adjoint) representation 
of the color group. In this way the Taylor expansion goes over covariant derivatives: 



[0,zMz) = J2 ¥ ki(D 22 ) k m- (2-54) 



k 



In what follows the [0, z]-factors are not shown for brevity, but they are always implied. 
Note that dropping the gauge links can be viewed as going over to the Fock-Schwinger 
gauge x a aA aa (x) = 0, A aa (0) = 0, or, alternatively, the light-cone gauge A 22 = 0. 

In addition, we have to replace ordinary derivatives by covariant ones in Eqs. (2.51). 
All relations which we have used to reduce an arbitrary combination of derivatives to this 
particular form hold true up to commutator terms [D aa , D"P] which can be expressed 
in terms of gluon field strength. Such terms contain two or more fundamental light-ray 
fields and do not affect the proof of completeness in the one-particle sector which we are 
considering at present. 

We still have to check, however, that the replacement d — > D does not spoil trans- 
formation properties of the basis fields (2.51). The special conformal transformation for 
the gauge field A takes the form 

i \^-aa.Ap^ (x) = X a -yXry^,d Anp(pc) -\- 2 {x a aApp(x) %flpA a Qt(x)^ 26^6^^ (l ( i) ) , 

(2.55) 

and for the components of interest becomes 

i[K 22 A ll }{x) =x 2 ^x j2 d^A li (x) , 

t[K 22 A 2i }(x) = {x^dr* + 2x 22 ) A 2i (x) - 2x 2i A 22 (x) , 

i[K 22 A 12 ](x) = (x 2 ^x l2 + 2x 22 ) A 12 (x) - 2x 12 A 22 (x) . (2.56) 

Using these expressions it is easy to check (by induction in m and j) that, e.g. for 
■0^! m (x) = (D 2 ^) 2 ' 7 2 (D 1 ^) 2m -?/' + (a;), the transformation is 

i [K 22 , ^ m ] W = (x 2 ^ 2 Sr> + 2jx 22 ) ^ m (x) + x 12 Q , (2.57) 

where Q is a light-ray operator containing some combination of the chiral field ip, gauge 
field A and derivatives. This inhomogeneous term vanishes on the light-cone, x = zn, so 
that one ends up with 

S + i>i m (z) = (z 2 d z + 2jz)^ m (z) , (2.58) 

as required. 
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The complete basis of one-particle light-ray operators for chiral quark and self-dual 
gluon fields in QCD contains seven fields: 



x+ 

X- 



++ 



./'- 



(1,771 

+ - 



{D^f-\D li ) 2m ^ 1 {z), 

(D 2i ) 2j 3 (D li ) 2m f 1 i(z) , 
(D^ 2j -\D^) 2m Uz), 
(D^ 2m f 12 (z). 



(2.59) 



The field carrying the superscript j transforms according to the representation T J of 
the SX(2,R) group, see Eq. (2.41). Note that ordering of the covariant derivatives in 
(2.59) does not affect the transformation properties. The twist E and helicity H take 
the following values: 





= (2j + 4mT l)^2' m) , 


7-1 _(i',m) 

Ex± 










= ± (23 - 1 =f I) ^i m) , 


TT -(ji m ) 

Hx± 









= (2 + 4m)/ 



(1,771) 



T 2 j - 1 T 



I 0; m ) 



=0 



(2.60) 



(2.61) 



The basis fields in the antichiral sector can be defined as ip 
for all other cases. 



(s,m) 



(s,m) ^ 



and similarly 



The proof that the fields in (2.59) form a complete basis in the one-particle sector 
essentially follows the above discussion of a chiral field. To this end one can consider 
the derivatives as commuting ones and assume that the fundamental fields satisfy "free" 
EOM . One has to demonstrate that all possible combinations of derivatives acting on the 
self-dual strength tensor can be reduced to the combinations appearing in (2.59). This 
can be achieved by inspection. The first step, as above, is to get rid of all pairs D 12 D 21 
replacing them by D 22 D 11 — > d z D u , and then check that all remaining combinations 
can be rewritten in the desired form, e.g. 



(D 12 ) 2k (D u y m f 12 (z) - -d z {D^ 2 Y K -\D^Y m f 22 {z) + 0{f 2 ) 
{D 2i ) 2k {D^) 2 ™ f 12 {z) - —{D 2i ) 2kl {D li ) 2m+1 f 22 (z) + 0(f 2 ) 



, 2m 



, 2Jfc— 1 . 



, 2m 



(2.62) 



etc. 
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Finally, taking a color-singlet product of the basis fields defined in Eq. (2.59), 
$i' m = {^ m , . . . , }, and their antichiral counterparts, 5>- ?,m , at different light-ray 

positions z%, . . . , z^, one obtains a complete basis of gauge- invariant iV-particle operators 

0(z u ...,z N ) = &^ mi { Zl ) . . . & N ' mN (z N ) (2.63) 

that transform according to the representation T- 71 ® . . . ® T^ N of the collinear conformal 
group SL(2, K) and serve as generating functions for towers of the local operators of twist 
E = Ei + . . . + En- If E > N then these operators get mixed under renormalization 
with the operators of the same twist E and the number of fields ranging from iV to E. 
Hence the mixing matrix has a block-triangular form. The anomalous dimensions are 
determined by the diagonal blocks only, the off-diagonal blocks are, however, important 
for the construction of multiplicatively renormalizable operators. The premium and main 



rationale for using the conformal basis (2.59) is that the SX(2,R) symmetry imposes 



severe constraints on the form of the kernels and also allows one to apply many of the 
technical tools that were developed earlier for quasipartonic operators. The explicit 
construction of this basis presents one of the main results of this paper. 

Since the maximum light-cone spin projection coincides, obviously, with the Lorentz 
spin, quasipartonic operators have definite geometric twist T = E = N. On the other 
hand, non-quasipartonic operators contain both T = E contributions and those with a 
lower twist, T < E. Operators with different values of T do not mix. Thus, introducing 
operators with different geometrical twist would bring the mixing matrix in the block- 
diagonal form at the cost, however, that the SL(2, M) symmetry of the kernels is lost. A 
better strategy is to separate the (highest) geometric twist of interest by imposing the 
appropriate symmetry conditions on the solutions of the renormalization group equation 



for the operators in (2.63) and maintain the SL(2, R) covariance. 

For illustration, let us consider a simple example: renormalization of twist-3 operators 
that one encounters in the study of chiral odd pion distribution amplitudes (32]. The 
complete set includes in this case three E = 3 light-ray operators 

O l (z 1 ,z 2 ) = x+(zi)ip-(z 2 ) , 

2 (z 1 ,z 2 ) = x~(zi)ip + (z 2 ) , 

3 { Zl , z 2 , z 3 ) = x+(z 1 )f ++ {z 2 )Mz 3 ) (2.64) 

that transform according to the representations T J=1 ®Ti =l l 2 , T J=1//2 (&Ti =l and T J=1 <g> 
yi=3/2 rpj=i Q f ^Yi e co ]}j near conformal group SL(2, M), respectively. The renormaliza- 
tion group equation can be written, schematically, as 

where the kernels H. 2 ^ have simple SL(2,M.) transformation properties. The operator 
03 is quasipartonic with T = E = 3, the other two are non-quasipartonic and do not 
have definite geometric twist. 
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For example, consider a generic local operator obtained by the Taylor expansion of 
the light-ray operator 0\ 



X a X c 



D. 



It is symmetric in all dotted indices, but does not have definite symmetry in the undotted: 
(3 can be either symmetrized or antisymmetrized with (either) one of the a, ati, . . . a^. 
These two possibilities correspond to picking up contributions of different geometric twist 
T = 2 and T = 3, respectively Going over to local operators is in fact not necessary 
as the separation of contributions of different symmetry can be achieved in the nonlocal 
form: 



^1,2 



It is convenient to introduce 



,«/3 







d 



■Oi, 



± (z 1 ,z 2 ) = 1 {z 1 ,z 2 )±0 2 (z 1 ,z 2 ) 



The operator 0_ can be written as 



(2.66) 



(2.67) 



(2.68) 



and is pure twist T 
after some algebra 



3, whereas for + inverting the relations in (2.66) one obtains 



+ (z u z 2 ) = /i a ^ j dr O t =\tz x , tz 2 ) - (/iA) J* dr r O t =\tz u rz 2 ) . (2.69) 

Using O- and (9^ =3 as basis fields instead of 0\ and 2 one can avoid the contamina- 
tion by twist-two operators altogether. The problem is that the renormalization group 



equation will in this form involve linear combinations of the kernels (2.65), of the type 
Tin ±7^12 etc., with different SL(2, R) transformation properties. 

In the present case the problem can be simplified drastically using the operator 
identities [32] that allow to rewrite both two-particle operators Q^ =3 and Q3 =3 in terms 
of O3 (up to a local term x a (0)V'a(0) [32]) so that they do not need to be considered 



separately. In this way the matrix renormalization group equation (2.65) is reduced to 



the single term 7^33. Unfortunately, a similar reduction to the quasipartonic sector does 
not hold in the general situation. 



3 Complete Operator Basis for twist-4 

After this general discussion we proceed to the systematic study of the renormalization 
of QCD baryon operators of twist-4. From the theory side, this is the simplest example 
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where non-quasipartonic operators enter nontrivially and cannot be excluded by using 
EOM; a generalization of the approach of Ref. [23J to such situations is our primary 
goal. The main application of twist-4 baryon operators to QCD phenomenology has 
been to the studies of hard exclusive reactions involving helicity flip, for example the 
Pauli electromagnetic form factor of the nucleon F 2 (Q 2 ) using pQCD factorization [31] 
or light-cone sum rules [35J. The necessary nonperturbatuve input is given in this case 
by the three nucleon (proton) distribution amplitudes $ 4 , and S 4 defined in Ref. [36] 
as matrix elements of twist-4 three-quark operators. They can be represented in spinor 
notation as follows: 

{^ k r+\z l ) X ^{z 2 )^ d l k {z,)\P) =- -^X)m N Nl J Vxe-^^'^x), 

(0\e^ k xf{z x ) i>y{z 2 ) ijf(z 3 ) \P) = - lfji\) m N Nl j Vx e"^ E^V^x) , 

{^^{z^iz^f fa)\P) =- ~(ji\)m N Nl J Vxe-^ n ^ x ^E 4 (x). (3.1) 

For comparison, the leading twist-3 distribution amplitude can be defined as 

(0\e^ k rA^)x U + j (^f(z 3 )\P)= \{pn)Ni j Vx e*** 2>** 8 (*) , (3.2) 

and only involves "plus" quark fields. In addition, we introduce three independent twist-4 
four-particle distribution amplitudes involving a gluon field: 

(0|z^7(^)x^(^)[/ ++ (^)^(^)] fc |P) =\m N {pnfNl j Vxe~^^ $f (x) , 

(0\ige^ k xf(zi) [f ++ {z,W + {z 2 )]^f{z,) \P) = -m N {pn) 2 Nl j Vxe~^ v^( x ) , 

(0\ige^ k [f + 4z^l(z^] i ^\z 2 )i>i k (z 3 )\P) = -m N {pn) 2 Ni j Vx e -'0»>£** E° (x) . 

(3.3) 

Here the first and the second superscripts of "0 s ' 1 {x a ' 1 ) are the flavor, a — u,d, and color, 
i = 1, 2, 3 indices of the quark field, respectively; m N and p^ are the nucleon mass and 
momentum, and = (1/2) (1 ± 75) N(p) is the antichiral (chiral) part of the Dirac 

spinor. In order to keep the auxiliary spinors X,fi dimensionless we choose 

n = A (g> A, n = fx <g> JJl . 

The distribution amplitudes depend on the set of parton momentum fractions x = 
{xi, . . . , x n } and the integration measure is defined as 

J Vx = J dxi . . . dx n 5(l — 'y~]xk) . (3.4) 
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We hope that using the same notation j T>x for the three-particle and the four-particle 
integration measure will not create confusion. We keep the factors (/xA) = 1 on r.h.s. of 



Eq. (3.1) to maintain the balance between the spinors /i and A on the both sides. 



The scale dependence of the distribution amplitudes is driven by the renormalization 



of the nonlocal light-ray operators on the l.h.s. of Eq. (3.1 ). Our first task is to construct 
the complete operator basis. Note that S4 involves chiral quarks only whereas $4 and \l/4 
involve both chiral and antichiral fields. Since chirality is conserved in QCD perturbation 
theory, there is no mixing and the two cases (pure chirality and mixed chirality operators) 
can be considered separately. 

3.1 Chiral operators 

The distribution amplitude S 4 is related to the matrix element of the operators with 
collinear twist E = 4 and helicity H = 1/2 

Q1O1, z 2 , z 3 ) =e<"ty*(* 1 ) i> b + J (z2) , 
Q 2 ( Zl , z 2 , 23) =e«Vrt*i) ^ j (z 2 ) r+\z 3 ) , 

Q 3 ( Zl , z 2 , z 3 ) =e^ k r+\zi) i> b + J (z 2 ) r-\z 3 ) . (3.5) 
For the discussion of the renormalization it is convenient to consider the three quark fields 



of different flavor in which case the three light-ray operators in (3.5) are independent. 
We do not assign flavor indices to Qi assuming that the flavors are always ordered as in 
the above expressions. 



The three-quark operators in (3.5) mix with quasipartonic operators containing an 



additional gluon field F + „^ = — (/xA)/ ++ 

G 1 (z h z 2 , z 3 , z A ) =igj*(ji\) if ++ (z*WM)Y t+fa) r+ k (z 3 ) , 
G 2 ( Zl ,z 2 ,z 3 ,z 4 ) ^ige^i^rAzt) [f ++ (z A )4j b + (z 2 )]^f (z 3 ) , 

G 3 (z h z 2 , z 3 , z A ) =tge^ k (fi\) ^ j (z 2 ) [f ++ (z,)r + (z 3 )} k , (3.6) 
which are, however, not all independent because of the identity 

Gi(zi, z 2 , z 3 , Zi ) + G 2 (zi,z 2 , z 3 , Z4) + G 3 (zi, z 2 , z 3 , Z4) — . (3-7) 

Recall that ?/>_, ip + and / ++ have conformal spins j = 1/2, j = 1 and j = 3/2, respec- 
tively. Hence Q\ transforms according to the representation <g> ® T^ 1 ', of the 
SL(2, K) group (and similar for Q 2 , Q3), whereas G { belong to the T^^T^^T^^T^^ 
representation.. 

For the matrix elements of the operators Qi and Gi between the vacuum and proton 
state there are more relations that follow from identity of the two u-quarks in the proton 
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and the requirement that the nucleon has isospin 1/2. Let 

q i {z 1 ,Z2,z 3 ) = (0\Qi{z h z 2 , z 3 )\P) , 
gi{zi, z 2 , z 3 , z 4 ) = (OlGiizt, z 2 , z 3 , z 4 )\P) , 

where we put a = b = u and c = d. One finds 

q 2 (z u z 2 ,z 3 ) = q x {z 2 ,z x ,z 3 ) , 

q 3 (z 2 , z 3 , zi) = -giOi, z 2 , z 3 ) - qi(z u z 3 , z 2 ) 



(3.8) 



(3.9) 



and similarly 



92\z\i z 2 , z 3 , z 4 ) — g\(z 2 , zi, z 3 , Zi) , 

93( z 2, z 3 , Zi, z 4 ) = —g^zi, z 2 , z 3 , Z4) — gi(zi, z 3 , z 2 , Z4J . 



(3.10) 



Taking into account the identity (3.7) it follows that the remaining independent function 



gi satisfies the symmetry relation 

g x {z u z 2 , z 3 , z 4 ) - gi{z 3 , z 2 , z-l, z±) = gt(z 3 , z u z 2 , z&) - gi{z 2 , z u z 3 , z 4 ) . (3.11) 
Going over to the momentum fraction representation, we define a new twist-4 four- 



particle proton distribution amplitude S4 in Eq. (3.3). 



3.2 Mixed chirality operators 

The distribution amplitudes $4 and ^4 are given by the matrix elements of the operators 
with collinear twist E = 4 and helicity H = —1/2. There exist three independent three- 
quark operators with these quantum numbers: 

Q 1 (z 1 , z 2 , z 3 ) =e^ k r-\zi) l^teO xf{z 3 ) , 
Q 2 ( Zl , z 2 , z 3 ) =e^ k i>?(zi) ^2) xf(z 3 ) , 
Qs(z 1 ,z 2 ,z 3 ) J- e ^^\ Zl )^{ Z2 ) [xTr\z 3 ) , 



(3.12) 



where x+ 2 = X+^ 2 '°' ) — —(l J 'DX)x+ = —D ^x+i c f- Eq. (2.59). In addition there are 



three operators containing an extra gluon field 

z 2 , z 3 , z 4 ) =ige^ k (fiX) [/ ++ (^ 4 )^(^i)] i V^fe.) xf(z 3 ) , 
<3 2 (z 1: z 2 , z 3 , z 4 ) =ige ljk (//A) ^+ l (z x ) [f ++ (zi)^ b + (z 2 )] 3 x+(z 3 ) , 
g 3 ( Zl ,z 2 ,z 3 ,z A ) =ige^ k ^X) ^1)^2) Ih+i^xK^T , (3.13) 
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which, again, are subjected to the constraint 

£1(21,22,23,24) + £2(21,22,23,24) + £3(21,22,23,24) = 0. (3.14) 
Thus there are two independent distribution amplitudes in this case, and ^4, which 



can be defined as in Eq. (3.3) 



The operators 21(22,23,2:1) and Q 2 (zi,z 3 ,z 2 ) with the choice of flavors b = d and 



a = c = u enter directly the definition of \I/4 and $4 in (3.1). For practical applications 
it can be convenient to introduce an additional distribution amplitude related to the 
operator Q 3 

^(0|e^r(2i)[xf P j (*2)V? (z 3 )\P) =\^\){pn)m N Nl j Vx e 4 ™ E ^ D 4 (x) , 

(3.15) 

The amplitude _D 4 is not independent and can be reduced to the combination of other 
amplitudes using EOM. For the matrix elements in coordinate space 

^(21,2:2,23) = (0\Q k {zi, z 2 , z 3 )\P) , ^(21,22,23,24) = (OlQkiz^ z 2 , z 3 , z 4 )\P) , (3.16) 

this relation reads 

/ s d / \ 9 . . 

Cp 3 {Zi, 2 2 , 2 3 ) = — <^i(2i, 2 2 , 2 3 ) + ^^2(21, 2 2 , 2 3 ) 

1 * 



, dr (213 (21, 22,2:3, z T l3 ) +Z23<P 9 2 (zi, 22,23,223)), (3.17) 

where we use the notation 

z ik = Zi-z k , f=l-r, z T ik = Zif + z k r. (3.18) 

Going to the momentum fractions one derives 

D 4 (x 1 ,x 2 , x 3 ) =x 3 $ 4 (a;i, x 2 , x 3 ) + xi^ 4 (x 2 , x lt x 3 ) 

1 / f X2 dx f X3 dx 

+ / — $i(x 1 ,x 2 - x,x 3 ,x) - / — ^l(x!,x 2 ,x 3 - x,x) 
^ \Jo x Jo x 

f Xl dx [ X3 dx \ 

+ / — ^{(x 2 ,xi - x,x 3 ,x) - / — $?l(x 2 ,xi,x 3 -x,x) ) . (3.19) 

Jo x Jo x J 

4 Renormalization Group Equations 



4.1 General properties 

The matrix of anomalous dimensions for a set of local operators Oj is defined as 

d_ 
dfj, 



llk = Z-^—Z ]k , [0]i = Z^O?, (4.1) 
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where [0]i, Of are renormalized and bare operators, respectively. The renormalized 
operators satisfy the matrix RG equation 



^ + ^)g- g +lik(9) ) \o\, 



o. 



(4.2) 



which, equivalently, can be cast in the form of an integro-differential equation for the 
generating functions (light-ray operators): 



Here 7 is an integral operator, which we write as 

a. . 



zn) = ■ 



7 



2vr 



e. 



(4.3) 



(4.4) 



where a s = g 2 j^n. In what follows we will refer to the evolution kernel EI as the 
Hamiltonian. It is determined by one-loop counterterms to the nonlocal operator. 

For the both cases of pure chirality and mixed chirality operators we will be dealing 



with three light-ray three-quark operators, (3.5) or (3.12), and three operators with an 



additional gluon field, (3.6) or (3.13), respectively. Thus EI is in both 6x6 matrix 

which can be written as 



. H, EI g5 
H " ( W 9 



(4.5) 



Each block, H 3 ,H 9 ,H 99 , is a 3 x 3 matrix where the entries are integral operators. The 
diagonal blocks, M q and M g , are given by the sum of two-particle operators, whereas the 
off-diagonal block M. qg describes 3^2 transitions. The explicit expressions will be given 
below. 

One of the three light-ray operators involving a gluon field, e.g. Gs(zi, Z2, z%, Z4) 



(3.6) and Gs(zi, z 2} z 3} Z4) (3.13), can be excluded from consideration with the help of the 



operator identity (3.7). Let 



f,chiral / 



(5) 







0,2(21, ■ ■ 




Qs(zi, ■ ■ 




Gi(zi, . . 


.,z 4 ) 


\G 2 (z u .. 


■ , z 4 )J 



(4.6) 



be the vector of the remaining five independent chiral operators, and similar for mixed 
chirality. This vector satisfies the RG equation with the modified Hamiltonian 







(4.7) 
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where EI is a 5 x 5 matrix such that 

[fi g ] lfe = \Bg\ik , i, k = 1, 2, 3 

[M g ] lk = [M g ] lk -[M g } l3 , i,k = l,2 

[M qg ] ik = [U qg ] ik -[W qg } i3 , i = 1,2,3, A; = 1,2 (4.8) 

The nonlocal operator O(z) can be expanded over a complete basis of local operators 

0(z) = Y,^N, q (z)0 N , q , (4.9) 

N,q 

where operator 0^, q has canonical dimension iV + 9/2 and q enumerates independent 
local operators of the same dimension. If a local operator Ojv, 9 satisfies the RG equation 
(fidn + P(g)d g + 7Ar,g) Ojv, g = then the corresponding "coefficient" function ^f NjQ (z) is 
the eigenvector of the Hamiltonian EI 

[EM N>q ](z) = E N ^ Ntq (z), (4.10) 

and ^N, q = (o: s /27i)EN jq so that 

where (3 = U/3N C - 2/3n f . 

The coefficient functions of local operators ^N, q {z) are homogeneous polynomials of 
three variables 

E ^ZA^4 3 , (4.i2) 

fcl,...,fe 3 
k 1 +k 2 +k 3 =N 

for the quark components, i = 1,2, 3, and four variables 

= E £ 4 2 4 3 4 4 (4.13) 

ki,...,ki 
fci+...+fc 4 =JV-2 

for the quark-gluon components, « = 4, 5. 

4.2 Conformally invariant evolution kernels 

The structure of the Hamiltonian is severely constrained by conformal symmetry. It 
follows from the group theory that a nontrivial operator mapping the representation 
T jl ® T J2 to T %1 <g> T* 2 only exists if the difference i\ + i 2 — ji — j 2 is an integer. If 
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^1+^2 = ji + 32 a conformally invariant operator K can be written in the form (see 
Ref. [37] for details) 



a/3 
a/3 



21 J ) 

(4.14) 



where the notation follows Eq. (3.18); the function k(x) is arbitrary. 

The two-particle one-loop kernels fall in two groups. The kernels in the first group 
involve a single integration, k(x) ~ 5(x), and their form is completely fixed (up to 
pref actor) by the conformal spins of the fields. For the case that the conformal spins are 
conserved, %i = j± and %i = j%, one important example i^j] 



1 da 
a 



a^-Xzi, z 2 ) - ip(z« 2} z 2 )\ + a^- 1 [ip{z u z 2 ) - ip(z Xl . 



(4.15) 



This structure is specific for gauge theories and arises from the diagrams involving the 
gluon field in the light-like Wilson lines (in covariant gauges). 

Another kernel of this type is the "exchange" kernel, 7i e . It maps T^®T^ 2 — > T^ 2 ®T^ 
and is fixed by conformal symmetry up to a prefactor [3?]. Assuming that j\ > 32- 



[HtM(zi, Z2) = / da a 2 * 2 ' 1 a 2 ^ 2 ^ 1 y(z^ z 2 ) . 



(4.16) 



The two-particle kernels (4.15), (4.16) depend implicitly on the conformal spins of the 
fields which they act on; this dependence will always be implied. 

The kernels (4.15) and (4.16) are both known from the studies of leading twist oper- 
ators. In addition, we will need 



[ni 2 ip](z u z 2 ) = v{z™ 2 ,z. 



21, 



daa 2h - l a 2j2 -\(z^z« 2 ) 



(4.17) 



which corresponds to the "diagonal" mapping without the spin exchange: T jl £g> T J2 — ► 
T jl ® T j2 with k(x) = 5(1 - x). 

The kernels in the second group retain both integrations as in (4.14) and usually 
involve a theta-function k{x) ~ 9(1 — x) or k{x) ~ 9(x — 1) which restricts the region of 
integration to a + (3 < 1 and 1 < a + f3, respectively. We define two more kernels 7i ± by 



[Htt<p](zi, Z2) 
[H^izt, z 2 ) 



da 



d(3a 2 ^ 2 B 2n - 2 ^Zi2,4i) 



da I d(3a 2]1 ~ 2 (3 2]2 ~ 2 ^(z^z p 21 ). 



(4.18) 



(4.19) 



^Here and below the subscripts Ttik, i, k = 1, 2, 3, 4, indicate that the kernel acts on the coordinates 
of the «-th and fc-th partons (particles) . 
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The off-diagonal block M. qg maps the SL(2, M.) representations with a different number 
of fields T 1 ® T 1 ® T 3 / 2 -> T jl ® T 3 ' a , where either j l = 1/2, j 2 = 1 or j x = 1/2, j 2 = 1 
so that in the both cases j% + j ' 2 = 3/2. The general form of the corresponding kernel 
consistent with conformal symmetry is 

[Rf](z u z 2 ) = 4 ^ da J" dp J 1 dr^-^-V fe ^ /(*&, 4, 4) • (4.20) 

where r(x, y) is an arbitrary function. In one-loop diagrams the following kernels appear: 

-1 r i 

dp 

a 



[ V i2(3) /] (*i> ^2) =4 / da j dp°^- f(z" 2 , z 2 , 4) 



pi pa 

[ v i ( 2( 3 ) /] (*i> *a) =4 / da / d/3 /3 /(z? 2 , z 2 , 4) 
Jo Jo 



[V 1 ? ( 3)/](^i^2)=4 / d/9 / dry ^-(^-2^) m,4,4) 



J/3 

[ v i2(3) f]( z i> z i) =z i2 dpp z 2 , 4) + ^ ^ ^7/(^1, 4> 4; J > 

/*! /»1 

\y$(3)f\( Z l> **) = Z 12 / rfa / d ^ / d lPf( Z l2i Z ln Z 2l)i 
JO Ja Jo 

[V<$ (3) f}(z h z 2 ) =4 /'da f d/3 rdjPf(z? 2 ,zl,z^), 
Jo Jo Jo 

[ v i ( 2(3)/](^i^2) =4 / d« / d/3 / d7/5/(^ 2 ,4,4) , 
Jo Jo Jo 

which correspond to the choices 

r (1 \x,y) = 8{x)9{x/y- 1), r {2) (x,y) = 5{y)9{\ - x/y), 

r®(x, y) = S(x) 6(y - 1) Q - 2V r^x, y) = 0(1 - y)<f(x/j/) (1 + %)), 

r (a) (x,?/) = 9(l-x)9{x/y-l), r (b) (x,y) = 6(1 - x)9(l - x/y), 

r®{x,y) = 9(l-y)9(l-x/y). 



4.3 The scalar product 



We will be looking for solutions of the Schrodinger equation (4.10 ) that are polynomials in 
light-cone variables 2%. The scalar product on this space can be constructed as follows. 
First of all, we allow z^ to take complex values. It is convenient at this point to go 
over from the SL(2,R) group to 577(1,1) which has the same algebra. In particular, 
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the SU(1, 1) generators have the same form. The SU(1, 1) invariant scalar product is 
defined as [38] 



D jZ f x {z)h{z), 



??>- 2 d 2 < 



\z\<\ 



71 



(4.32) 



where j is the conformal spin, j > 1/2, and functions fk(z) are polynomials in z; f(z) = 
(f(z))* stands for complex conjugation. For the special case j = 1/2 the SU(1, 1) scalar 
product takes the form 



</i,/a> 



J=l/2 



1 

2^ 



2tt 



d<p /i(e*) / 2 (e^) . 



(4.33) 



One easily finds ||^ n || = (z n ,z n ) 



r(2j)n! 

— r-. For polynomials of several variables 

r(2j + n) 

f(zi, . . . , z n ) which belong to the tensor product T- 71 ® . . 
given by 



T^ n , the scalar product is 



(/i,/a> 



n / D i*^* -fr^ 1 ' 

v fc = Vi**i<i / 



^n) fl{. z lf--t z n) 



(4.34) 



It can be shown that the diagonal blocks (three-particle and four-particle) of the evolution 
Hamiltonians are self-adjoint with respect to the following scalar product: 



(* 1 |* 2 > = ^(*f ) |n tt *S* ) > 



(4.35) 



where the matrix Q is 



n 



fa x 0\ 
a 2 
a 3 
2 1 
\0 1 2/ 



(4.36) 



The coefficients aj = a for the chiral operators and a>i = a 2 = 2a 3 = a for the mixed 
ones, where a > is an arbitr ary co nstant. We put a = 1 in what follows. I t is t acitly 
implied that (^Ifi^^), in ( 4.36 ) is given by the SU(1, 1) scalar product J4.34) with 
the spins equal to the conformal spins of 

We would like to emphasize that the choice of the scalar product depends on the 
problem under consideration, see, e.g. [39J for a discussion. 



4.4 Chiral operators 

We are now in a position to write the explicit expression for the Hamiltonian that 
governs the RG equation for the chiral operators (3.5), (3.6). The quark block can 
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be written in terms of the two kernels denned in (4.15) and (4.16): 

1 



1 + 



(4.37) 



with 



n 




(4.38) 



where 



TT %JV I n_/v I n/v 



23 



L 31 



(4.39) 



As discussed in Sect. 3 the Hamiltonian W q describes renormalization of the three quark 
operators of both geometrical twist T = 3 and T = 4. On the other hand, the Hamil- 
tonian for the twist-3 operator is known from Ref. [8], where it is called H 3 / 2 . 
Thus one obtains a nontrivial consistency condition^] 



(H^) fcl + (HfW% 2 + (Hj^)fc3 = H 3/2 



k = 1,2,3 



(4.40) 



which, alternatively, can be used to restore M q from the known result [8] for H 3 / 2 . 

The "gluon" block M. g describes the renormalization of four-particle quasi-partonic 
operators and can be restored from the results existing in the literature [HI [231 40J. 
Though the baryon operators in question make sense for N c = 3 only, it is convenient to 
separate the terms with the different color factor: 



ETfW = N c + e( 0) + — 



21 

y 



(4.41) 



where the last term (a constant) stands for the self-energy type contributions. For the 
diagonal elements we find 



-!)■ 



kk 



kk 



kk 



_1_/ v q 



fc4 ' 



z 7~Lk+l,k-l + ^fc+1,4 + ^Gc-1,4 ^ \ H~k+l,4 + ^fc-1,4 + ^fc+1,4 + ^Jfc-1,4) 



^12 + + ^Sl — ^ At 



L 23 



31 



L fc4 



(4.42) 



In these expressions k = 1, 2, 3 and k, k ± 1 appearing in the subscripts of the kernels 
on the r.h.s. refer to arguments of the quark fields. If/c + l = 4orA; — 1 = the 
corresponding subscript should be put to 1 or 3, respectively. 



H3/2 on the r.h.s. of (4.401 is given by the sum of three 'V-type kernels for the conformal spins 
pairs {ji,jk} = {lj 1} whereas H° on the l.h.s. contains one kernel with {ji,jk} — {1, 1} and two kernels 
with {ji,jk} = {1,1/2}- The "e"-type kernels appearing in the non-diagonal entries (H g )i 2 + (H g )i 3 
correct for this difference. 
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The off-diagonal matrix elements i, k = 1, 2, 3, i ^ k, read 

[Hf] 



ifc — n ik n k4 z rl fc4 

-2HI 



1 

2 ' 



The identity in Eq. (|3.7|) implies the relations 

[H 9 ] 13 + [H s ] 21 - [M g } 23 = [M g ] 33 - [H s ] 31 , 

[h 9 ] 23 + [H s ] 12 - [H fl ] 13 = [e 9 ] 33 - [e 9 ] 32 . 



[H 9 ]n 

ft] 22 



Finally, for the off-diagonal quark-gluon block we find 

1 



ET 



chiral 
99 



1-/ 



chiral 



with 



[TV 



chirali 

99 J 



A'/r 



1 'v (1) +v (1) 



chirali 

qg 



\ik 



N, 



V (1) + V 



V. 



(2) 

fc,fc+l,(4) 



V. 



(2) 

fc,fe-l,(4) 



(2) 
ifc(4) 



where, as above, the subscripts k,k ± 1 take the values 1, 2, 3. 



4.5 Mixed chirality operators 



For the mixed chirality operators (|3.12|), (|3. 13 ) the quark block is 

1 



i j z_ 1 o/wx 



where 



/H + HL-H 



13 

e 
21 



L 2.3 



7i e 



H + — 7^ 



12 
23 ~ 



13 



^13^13 
z 23^ 23 



\zT 3 1 (l-2Hf 3 ) z£(I-2H d 23 ) R-2(nf 3 +nt 3 )+Sj 



Similar to Eq. (4.40) one can derive the consistency relations 

\a^] 11 + \E^] 12 -\E^] 13 d za 
[n^ x ] 33 9 23 -[e^] 31 -[e^] 32 



H 1/2 



dz 3 H l/2 , 



(4.43) 



(4.44) 



(4.45) 



(4.46) 



(4.47) 



(4.48) 



(4.49) 



where the Hamiltonian Hi/ 2 describes the renormalization of the leading twist-three 
nucleon distribution amplitude [5]. 
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The gluon block HT^ X ^ can be expanded in powers of N c in the same way as in (4.41 ), 



with 



kk 



--H v u -2(l-5 k3 )H 



kA > 



(4.50) 



CO- 



^23 + ^4 + n-zA - n 23 - zn 9A - zn 9A + -r 34 n 



L 24 



L 31 



L 23 



21 



L 24 



34 H-43 ) 



[w pm, (o)] 33 =w?2 + ^ + ^ _ 2 (h+ + nt 4 + + n M ) , 

i 



pfjV^x/,(o)] 



12 



H\2 T~^OA 27^94 



L 21 



21 



21. 



'1-21 ''•14 i" •^'I-iA 



-14 



11 



L i3 



[ H ^x/"(o)] 3j 



njv n_/v 0/+ i oi7 
- ™ j4 - ^j3 + ln 



1 

./'4 ~~ 2 ' 
(4.51) 



[Hf x/ ' ( " 1) ] fcfc = H v 12 + H\ 3 + H" 23 - H+ 3 - H+ 3 - 2(1 - 5 K3 )H M + hzP^Hlz , 



[HjM(-i) ]l2 = _ 2H - 4) 

[e ^X/-(-l)]. 3 = p 34H e 3; 



,(-1)1 



[ H ^x/,(-i)] 3 . 



-2H\ 



-2H 



14 5 



(4.52) 



where P 34 is the permutation operator Puf( z i, z 2 , z 3 , z 4 ) = /(z i, -22, Z a_, z 3 ), k = 1,2,3 
and j = 1,2. The entries [M^ xf ] ik satisfy the same constraint 14.441) . Finally, for the 
off-diagonal quark-gluon block we find 

1 



mixed 
19 



mixed 



where 



L ig\ jk 



2**9 



J,k = 1,2 



(4.53) 



\n-i 

\T-/ 



mixed] 
99 

mixed 
19 



?,k 



33 



V, 



(6) 



(a) 



-v, 



(3) 



"<2fc3 V (4) 3 (4) 2 fc3 ( 4 ) 2 fc3 ( 4 ) 

2 



4 



fc = 1,2 
1 



w +_ ± v (c) +-v (3) +-v (4) 

j3(4) + 3 ^3(4) + g ^3(4) + 2 ^3(4) 



(4.54) 



^chirai ig gj ven j n Eq. (|4.46|) and 



<-W 



[AW M ]i2 


= [AW W ] 


21 =0, 




-V a) 

~3 K i3(4) 




[ A ^<Ji3 


-V (a) - 

~ ^3(4) 


_V b) +-v (c) 

3 K i3(4) + 3 ^3(4) 



(4.55) 



for j = 1, 2. 
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5 Renormalization Group Equations II: Solutions 



The Hamiltonian (4.7) has block-triangular structure. In order to find its eigenvalues it 
is, therefore, sufficient to consider the diagonal quark and quark-gluon blocks separately. 



5.1 Chiral quark operators 
5.1.1 Permutation symmetry 



The three-quark Hamiltonian T-^ww j s given by Eq. (4.38 ) with the kernels TC V , TC e defined 
in (4.15), (4.16), respectively. It is easy to check that 7 f(v^ commutes with the generator 
of cyclic permutations 

[P,-H^]=0, V = P a ®P z , P 3 = l, (5.1) 

where P a permutes the quark quantum numbers and P z the quark coordinates, respec- 
tively: 



P« *2,Z 3 ) 
P z Z 2 , Z 3 ) = VWJzs, Zl , Z 2 ) . 



^\z 1: z 2 ,z 3 ) 
(0 



(5.2) 



The eigenfunctions of can always be chosen to have definite parity with respect to 

the cyclic permutations: 



£6{l,e a / 3 1 e" &/3 }. 



(5.3) 



The (vector) eigenfunction ^f^ a can be written in terms of the single function ifj$ a as 



^%n{ z l,Z 2 ,Z 3 ) 



\ £2 ^N,g( Z 3,Zl,Z 2 ) / 



(5.4) 



The eigenfunctions of different parity are orthogonal with respect to the quark part of 
the scalar product (4.35), (\l/ £ |\l/ e ') ~ 5 ££ /, whereas for the eigenfunctions of the same 
parity one gets 



(5.5) 



The scalar product on the r.h.s. is given by Eq. (4.34) for the spins j\ = 1/2, j 2 = j 3 = 1. 
In addition to the symmetry under cyclic permutations, the Hamiltonian %pH> com- 

(12) (12) 

mutes with the permutation operator V% 2 = Pa <E> Pz defined as 

P^¥'\z) = ¥ 2 \z) , P^¥ 2 \z) = ¥'\z) , P^V(zi,z 2 ,z 3 ) = ^{z 2) z u z 3 ) . 

(5.6) 
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Since Vi 2 V = V^Vw one concludes that V12 ^% q ~ q • Thus the eigenvalues of the 
Hamiltonian 7-^^ in the sectors with e = e ±l2n ^ 3 coincide. 

Furthermore, if ty e N is the eigenfunction of the Hamiltonian 7-^^ with the eigenvalue 
£ft , then ip £ N is the eigenfunction with the same eigenvalue of the effective Hamiltonian 

n(e). 

^N^U Z 2, Z s) = ^N,q^N,q( z l, z 2,Z3) , (5.7) 

where 

H(e) = H 3/2 - HUl - eP- 1 ) - H e 13 (l - e~ x P z ) . (5.8) 

Here 7^3/2 is the evolution Hamiltonian for the twist-3 chiral quark operator [8], 

3 

7^3/2 = U\ 2 + 7^23 + ^31 + 2 ' ( 5 - 9 ) 



where the three kernels TC!- k are all given by (4.15) for the conformal spins ji = jk = 1. 

As was explained earlier the Hamiltonian 7^w ; a li as TC(e), describes the evolution 
of chiral operators of both geometrical twist-4 and twist-3. The twist-3 eigenfunctions 
^^^(zi, Z2, Z3) which belong to the sector with P-parity e, have the same parity with 
respect to cyclic permutations of the coordinates P z alone, P 2 \E^~ 3 = e^^' 3 . From 



the explicit expression in (5.8) it follows, obviously, that on such functions Ti{e) reduces 
to 7^3/2, as expected. 

It follows that the eigenvalues (and eigenfunctions) of 7^^ can be separated in three 
symmetry classes corresponding to the renormalization group equation for: 

• twist-3 operators; 

• twist-4 operators with P-parity e — 1; 

• twist-4 operators with P-parity e = e ±t2n ^ 3 . 

In the application to the nucleon distribution amplitudes we are interested in the eigen- 



functions that satisfy the relations in (3.9) (otherwise, the matrix elements of the cor- 
responding local operators between the vacuum and nucleon state vanish). One can 
easily verify that the following combination of the eigenfunctions with e = e +2lw ^ 3 and 
e = e~ 2t7T / 3 has the required symmetry property: 

+ (5-10) 



The twist-4 eigenfunctions with e = 1 do not satisfy the relation (3.9). They are relevant 



e.g. for the distribution amplitudes of baryons with isospin 7 = 3/2. 
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5.1.2 Complete Integrability 



The Hamiltonian (4.38) possesses a hidden integral of motion. Let 

S ik = d k (z k - Zi) = (d/dz k )(z k - Zi) . 



(5.11) 



It is easy to see that S ik acts as the intertwining operator between the representations 

T h=l/2 T ji=l and T j k =l T j t =l/2. 



S ik T jk=1/2 <g> T k=1 = T jk=1 <g> T u=l/2 S ik . 



(5.12) 



We define two-particle (3x3 matrix) operators [Qf k ] % ' k ', where i < and z', fc' = 1, 2, 3, 
by 



[Q 



[Q 



ik\ ~ ^ki 



(5.13) 



and all other off-diagonal matrix elements being zero. For the diagonal matrix elements 
we put [Q%f = [Qf k ] kk = i and 



2 

[QW = \ + s lk 



[QikV 1 — n + Ski, 



for j different from i and k. For example, explicit expressions for Q 12 are 

S 12 



Q+=h+\ $, s 




Qr 2 = ^i+ |^2i o o 

1 5- 



21/ 



The two-particle Casimir operators Jf k can be written in terms of Q ik as- 

Jik = 2 Wife' ^jfc) — ^ 
The Hamiltonian Jj^w^ can fog represented in the form 



(5.14) 



(5.15) 



(5.16) 



n^ = n 12 + n 23 + n 31 + -, 



where 



n 



12 



'H\ 2 - 1 H\ 2 
^12 ^12 _ 1 

o o n\ 2) 



and similarly for H 23 ,H 3 i. 



(5.17) 



(5.18) 



§It is easy to see that the anticommutator jf k of the Qf k operators is SL(2,M) invariant whereas the 
Qf k themselves are not. 
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Figure 1: The spectrum of the conserved charge, Qz/(N 
chiral quark operators. 



3) 3 , for twist-3 and twist-4 



Finally, lefl 



Qs 



-\J> J 2 } 

2 V J l2i J 23\ 



The operator Q 3 commutes with the Hamiltonian 7^/"^: 



(5.19) 



[Q 3 ,H. 



q J 







(5.20) 



and defines, therefore, a nontrivial integral of motion. To prove (5.20) it is sufficient to 
show that 



[{n ik ),Q*) = i{Ji j -% 



(5.21) 



where j ^ i,k, and the pair-wise Hamiltonians (Ttik) are given in Eq. (5.18). This 
can be done by calculating both sides of the relation (5.21) in the conformal basis, see 
Refs. US] for the details. 

It is straightforward to check that the operator Q 3 commutes with the operator of 
cyclic permutations V and anticommutes with Vu, 



[Q3,V} = {Q 3 ,V 12 } = 0. 



(5.22) 



Eq. (5.22) together with (5.20) imply that all eigenstates of the Hamiltonian are double 
degenerate, except for the states which are annihilated by Qs i.e. Q-^ = 0, cf. |9J. The 
spectrum of Q3 is shown and compared with the corresponding spectrum of the twist-3 
conserved charge [H] in Fig. [TJ 



1 Note that [J? 2 , Jf 3 ] = [J*,, J*] = [Jf 3 , J? 3 ] 
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It was shown in Ref. [27] that the spectrum of one-loop anomalous dimensions of 
(anti)chiral composite operators in QCD, i.e. operators constructed from the (anti)chiral 
fields and their derivatives, coincides with the spectrum of a certain (integrable) SU (2, 2)- 
invariant spin chain. The Hamiltonian (5.17) can be viewed as the restriction of the 
general SU(2,2) spin chain Hamiltonian on the subspace with twist E = 4 and helicity 
H = 1/2. 



5.1.3 The spectrum of anomalous dimensions 



A short-distance expansion of the nonlocal operator O(z), see Eq. (4.9), runs over a 
complete set of local operators On,q including operators with total derivatives. It is clear 
that in order to find the anomalous dimensions the operators with total derivatives can be 
omitted. The operators without total derivatives can be singled out by their properties 



under conformal transformations: they transform according to (2.41) and are usually 



referred to as conformal operators. The coefficient functions ^N,q(z) corresponding to 
the conformal operators satisfy the following constraints 



(0L + =0. 



(2) 



v(3) 







(5.23) 
(5.24) 



that follows from the requirement that such operators correspond to highest weights of 
the corresponding representation. To get the first equation, we apply the operator P 22 



which generates shifts along the "plus" light-cone direction to Eq. (|4.9|). One gets 

(z)d. 



X)[(0i + 3, + 3 3 )*jv, 9 (2) 

N,q 



> J N,q 



N,q 



v N,q ■ 



(5.25) 



The r.h.s. of this identity contains only operators with total derivatives, hence the 
coefficients of the confo rmal operators on the l.h.s. must vanish. 

To derive Eq. (5.24) we apply the transverse derivative P 12 to the nonlocal operator 

Taking into account that 



of leading twist-3: 
z[P 12 ?/> + ](z) = —2d z ip_(z), one obtains 

3 



z 2 )^+ k (z 3 ) 



k=l 



(1_ 

dz k 





N,q, 



(5.26) 



where Qu{z) are defined in Eq. (3.5). Again, since the l.h.s. only contains operators with 



the total derivatives, the coefficients of conformal operators on the r.h.s. have to vanish. 

The spectrum of the Hamiltonian %pH> can be studied using powerful Quantum 
Inverse Scattering Methods (QISM) [UJ which is however beyond the scope of present 
paper. For this work we adopted a "brute-force" method, calculating T^w^ in the basis 
of functions 



e N , k (zi, z 2 ,z 3 ) 



k\(N-k)\) 



(5.27) 
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Figure 2: The spectrum of the Hamiltonian Bf M> = (1 + l/N c )H^. 
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N 



(5.21 



k'=0 



and diagonalizing the resulting (N + 1) x (N + 1) matrix (Ti^^)k'k numerically. The 



translation invariance of ejv.fc guarantees that the constraint in Eq. (5.23) is satisfied 
identically. 

The results are presented in Fig. [2j The spectra of twist-3 and twist-4 operators are 
shown in the upper left and upper right panels, respectively. The two lower panels show 
the twist-4 spectra in the sectors with £ — I and e = e l2n ^ 3 separately. In addition, the 
numerical values of the eigenvalues for N < 6 are collected in Table [2j 
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En,o 
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En,2 


En,3 
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4 * 
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3 
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7.724361* 



Table 2: Anomalous dimensions of twist-4 chiral-quark operators in units of a s /(2ir); N is 
the total number of covariant derivatives. The entries marked with an asterisk correspond 
to the operators with 'P-parity e = e±* 27T / 3 and the remaining ones to e = 1. All anomalous 
dimensions except for the lowest ones, En,o, for odd N = 2k + 1, are double degenerate. 



The general features of the twist-3 and twist-4 spectra are similar. For both cases 
all eigenvalues are double degenerate (see above), except for those corresponding to zero 
eigenvalue of the conserved charge Q3^N, q = Q3^N, q , Qz = 0. These special eigenvalues 
turn out to be the lowest ones in the spectrum and can be found explicitly. There are 
two series of such states, one in the twist-3 sector and one in the twist-4 sector. 

The twist-3 eigenstates with Q 3 = were studied in Ref. 0. They exist for even N, 
are invariant under cyclic permutations and have the energy 

e$-Jo(n) = ( 1 + jr) { 4 ^ (iV + 3) " ^ (2) ] " 1} ' N - even ' (5 - 29) 

where ip(x) is the Euler ^-function. The twist-4 eigenstates with Q3 = exist for odd 
N. They are also invariant and cyclic permutations (e = 1) and have the energy 



EZ=o( N ) 



1-^(2) 



N - odd. (5.30) 
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Figure 3: The spectrum of the Hamiltonian W^x = (1 + 1/N C ) H^ % . 



The appearance of (N + 3)/2 as argument of the Euler -0-function is characteristic for 
systems involving half-integer conformal spins, cf. [T2"l |4"2"] 

As is seen from Fig. [2] the upper parts of the twist-3 and twist-4 spectra are very 
similar to each other and are in fact interlacing. In particular the line of the largest 
eigenvalues is the same for both twists and is given by [9] 

£ max (A0= ( 1 + ^r) {61niV-31n3 + 67E + ^ + 0(l/iV)|. (5.31) 

The distance between the neighboring eigenvalues is in both cases 0(1/ N) in the upper 
part and 0(1/ In 2 N) in the lower part of the spectra, respectively. 



5.2 Mixed chirality quark operators 

The analysis of operators of mixed chirality goes along the same lines. The eigenfunctions 



of conformal twist-4 operators have to satisfy the constraint in Eq. (5.23) and in addition 



d r 



,(!) 







,(2) 



(5.32) 



instead of (5.24). The spectra of the operators of geometric twist-3 and twist-4 are shown 



in Fig. [3j They are non-degenerate in both cases, in difference to the chiral operators. 
The numerical values of the twist-4 eigenvalues for iV < 4 are collected in Table [3] 

The eigenvalues in the upper part of the spectra for twist-3 and twist-4 are very close 
to each other and to the corresponding eigenvalues for chiral operators, cf. Fig. [2] In 
particular, to the 0(1/ N) accuracy the line of the largest eigenvalues does not depend 



on chirality. It is the same for both twists and is given by Eq. (5.31). 



37 



N 







-2 


1 


2/9 


2 


2(14- V43) 




9 


3 


197-V5089 




45 


4 


3.706620 



En,i 
2 

32/9 



49-V73 
9 

589-^11161 
90 



JV, 2 



2(14+^43) 
9 

49+^73 
9 

6.634936 



N,3 



E 



NA 



1974 



'5089 



45 

589+^11161 
90 



7.858442 



Table 3: Anomalous dimensions of twist-4 quark operators of mixed chirality in units of 
a s /(27r); iV is the total number of covariant derivatives. 



5.3 Quark-gluon operators 

Anomalous dimensions of four-particle twist-4 quark-gluon operators correspond to the 
eigenvalues of the Hamiltonians Wp^f an d HWw f or the pure and mixed quark chirality 

cases, respectively. 

The Hamiltonian is given explicitly in Eqs. rt4.4lh-rt4.43b. It has the same 



symmetries as and its eigenfunctions can be classified by parity with respect to the 
cyclic permutation V of the three quarks (5.1), e = l,e ±l2n ^ 3 Due to the identity (3.7) 



we are interested in eigenfunctions which satisfy the restriction 



(5.33) 



The eigenfunctions belonging to the sectors with e 



e ±i2n/3 h ave the same eigenvalues, 



E 



E 



and are related to each other as ^% q { 



Pi 2^ 



l2*N,q\ 



As in the 



case of quark operators, in the applications to nucleon distribution amplitudes we are 



interested in the eigenfunctions of particular symmetry, cf. (|3.10|). One can easily verify 
that the combination ^ N ^(z) = (1 + P 12 ) 
properties. In turn, the eigenvalues with e 



[z) with e = e i2n / 3 has all the necessary 
1 are not degenerate and the corresponding 
operators (eigenfunctions) are relevant e.g. for the A-baryon. 

In order to calculate the spectrum we have used the basis of functions 



^N,k,m{Zl, Z21 Z3, Z4, 



,2a - Z4) [Z 2 



Z4) m (z 3 



z 4 



>N-2-k-m 



k\m\ (JV-2-Jfe-m)! 



(5.34) 



and diagonalized the resulting N(N — 1) x N(N — 1) matrix numerically. The translation 
invariance of tN,k.m corresponds to the restriction to conformal operators, cf. Eq. (5.23). 



The combined spectrum of three-quark and three-quark-gluon chiral operators is 
shown in Fig. [Zj (all parities) and in Fig. [5] ( e = e* 27r//3 only). Note that the quark-gluon 
spectrum is much more dense as for given A" there are O(N) quark and 0(N 2 ) quark- 
gluon eigenstates. Also, it is seen that for large A" the spectra overlap significantly. It is 



38 



35 

30 

25 

_ 20 
tr 

w 

10 
5 




I 



o 8 



8 8 



+ + 

+ 



a + 

+ + 

+ + 
+ 



+ 

+ + 



+ 

+ + 
+ + 



10 
N 



15 



20 



Figure 4: The combined spectrum of the Hamiltonians (crosses) and HJ^* (open 

circles), all parities. 



easy to show that for iV — > oo the eigenvalues lie within the bands 

1 R 1 R 

=fhiN<Efff <8lnN, ^-\nN < E^ f < 14 In AT . (5.35) 

Note that the lowest four-particle quark-gluon eigenvalue has the same logarithmic 
asymptotic as the lowest three-quark one, so that they are separated at most by a 
constant. This suggests that for most of the quark eigenstates there is strong mixing 
with the quark-gluon ones, see the next Section. 

The combined spectrum of three-quark and three-quark-gluon operators with mixed 
chirality is shown in Fig. [6j It is similar to the chiral case, except that all eigenvalues 
are nondegenerate. 



5.4 Multiplicatively renormalizable operators 

For the construction of multiplicatively renormalizable operators one has to take into 
account the off-diagonal quark-gluon blocks H^ iral and EI™ lxed for the chiral and mixed 
cases, respectively. Note that the complete (reduced) evolution Hamiltonian EI is not 
Hermitian, therefore its eigenfunctions are not mutually orthogonal. 

In order to find the multiplicatively renormalizable operators we adopt the following 
procedure. The nonlocal operator O(z) can be expanded over a complete basis of mul- 
tiplicatively renormalizable local operators (4.9) with the coefficient functions ^^ s {z) 
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Figure 5: The same as in Fig. |4j but for e = e t2n ^ 3 only. All eigenvalues are double- 
degenerate. 



that are eigenfunctions of the Hamiltonian EI 

cf. (4.10). Let $L be the eigenfunctions of the adjoint operator Hit 

[&^ N J(z) = E N ^ q (z). 
The functions $?N, q (z) and form a bi-orthogonal system: 



(5.36) 



(5.37) 



(5.38) 



Eigenstates corresponding to different eigenvalues are orthogonal with respect to the 
scalar product (4.35). The multiplicatively renormalizable local operators 0^, q can be 



obtained as the scalar product of the nonlocal operator with the eigenfunction of the 
adjoint operator 



v N,q 



(5.39) 



With the increasing N the expressions rapidly become very cumbersome so that in this 
work we present explicit results for N = 0, 1, 2 only. 
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Figure 6: The combined spectrum of the Hamiltonians H^ x (crosses) and H^ x ^ (open 
circles). 



To begin with, consider chiral operators. Let 

QVnteto) =e ii*[(„ . D) kl ^ a _Y [(n ■ L>) fe2 Vy j [(n • D) ki ip 
Qfateto) =e w*[( n . D)*^-]* [( n • D)*^-] 3 ' [( n . D) k3 ijj c ]k 



c ]fc 



+J ' 



c 1 k 



(5.40) 



and 



ig^ k (fiX) [(n • !>)*/++(»» ■ Df 1 ^? [(" ■ [(« ■ -D) fe3 ^+] 



\k-ijx ]k 



G (k u k 2 MM) =ige ijk^ [(„ . £))*i^aji . D) ki f ++ (n ■ D) k2 iP b + ] j [(n ■ D) k ^ c + ] k , (5.41) 

As discussed above, the eigenfunctions alias multiplicatively renormalizable operators 
can always be chosen in this case to have definite parity with respect to the cyclic per- 



mutations (5.3). We remind that the spectrum for e = e ±l2w ^ 3 is double degenerate. 



Any linear combination of the corresponding eigenfunctions satisfies the evolution equa- 
tion. We present the results corresponding to the operators of definite symmetry under 
permutation of the first and the second quark 



N,q 



[l±P 12 )n, q (5-42) 
which are more convenient for applications than those with definite e = e ±l2w ^ 3 . One 
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obtains 



chvral, + =g (000) + g (000) _ 2g (000) 

oi h r l > + =q{ 100) - -Qi 010) + gi 001) - -Q { 2 00) + Qf w) + QT l) + -Q§° 0) 

2 2 2 

o?r l ' + =Q ( i 00) + 2gS 020) + gS 002) - 4Qf 10) + 2qS 101) - 4Qi 011) + 2g?° 0) 



+ + ; - 4Q 2 7 - 4g^ u ; + 2Q V i> - 3Qr ' - 

7 ^(0000) 7_ 

12 Gl 12 



- 2gf 2) + 8gf 0) + 2Q< 101 > + 2gf 1} - ^ G ; oooo) - -g$° 00) 



*g,chiral 
V 2,0 



,+ = 3 ^(0000) + ^(0000) j ^ ^c^a/ = lg/ g ^ 



and 



NchiroJ,- _ n (000) _ n (000) 
V 0,0 ~Vl V2 ) 

jcMrol,- =3g (100) + l g (010) _ 2g (001) _ 1 Q (100) _ 3 g (010) + ^ g (001) 

2 2 



=gS 200) + \Q { r ] + ^gS 002) - 2 gS 101) - 4gS 011) - \qT 0) - Qf 0) - \qT 2) 



+ 4Q { 2 101) + 2gf J > + Jgf °) - Igf °) - 2 gf 01 > + 2gf J > 

7 ^,(0000) 7 (0000) 
~ 36 1 36 2 

QftcMroi - = 1 ^(0000) _ G (0000) j ^ £f£ Wra ' = jg/g _ (544) 

The multiplicatively renormalizable operators of the lowest dimension in the e — 1 sector 
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arc 



Ochiral,l 
1,0 

^chiral,la 



*chiral,lb 
V 2,l 



=g (010) _ g (001) _ g (100) + g (001) + g (100) _ g (010) 



Qr> _ g (002) _ 6g (110) + 6g (101) _ g (200) 



j (020) 

+ g 3 

=g (200) 



6Q 2 



(200) _ n (020) 
V3 



(020) 



(002) 



Qi 110) -Qr ij + 2Qi' 



(101) 



2 

(020) _ 1 

2 



Qr ) -Qr ) +2gr ) -Q 2 



The operators 



jjchirai,la 



1 

2 
and 



g (020) +g (002) + 2 g (110)_ g (101)_ g3 



(Oil) 

(Oil) 
2 

(Oil) 



(5.45) 



<schiral,lb 
"2,1 



have the same anomalous dimension Ef^ = 4, 



cf. Table [2j so that any linear combination of them is multiplicatively renormalizable as 
well. Note also that for N = 2 there is no quark-gluon operator, the reason being that 
for e = 1 the cyclic permutation symmetry picks up the combination which is forbidden 
because of the condition in Eq. (3.7). 

Next we consider the operators involving quark fields of mixed chirality. Let 

Q (kuk 2 ,k 3 ) = £jk ^ n . Jjf^aji ^ n . £>)fc2^]J . D) k3 X C +) k , 
Q (kxMM) =e ijk ^ n . D ^ay ^ n . £>)fc 2 ^& y . D) k:i X C +} k , 



(fcl,fc 2 ,fc3) 



[(n • £>) fcl ^P [(n • £>) fe2 ^] j [(« • i?)**^ 



(5.46) 



where x+ 2 = xV 3/2 ' 0) = -(jiD\)x+ = -D xX+, cf- Eq. (2.59), and 



g ( kl MMM) =ige ijk (/iA) [(n . D) k ir+] i [{n . D f,] ++{n . D f^ + y [(n . . (5.47) 

In the mixed case there is no additional symmetry and all eigenvalues are non-degenerate. 
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The multiplicatively renormalizable operators of the lowest dimension read 



mixed 
0,0 


= Q (000) 


- Q { r ] , 








mixed 
1,0 


= Q (100) 


+ Qi 010) - 


3 (ooi) 
2 1 


- qS° 0) - Qf 0) + 


3 n (ooi) 
2 ^2 , 


mixed 


= Q (100) 


- Q { r ] + 


1 n (ooi) 


- Q? 0Q) + Qf 0) + 


1 Q (001) + Q (000) ; 


mixed 


= Q (200) 


_ Q (020) _ 


2 (002) 

3 1 


- 2ei ul) + 3ef u) 


- q? 00) + qT 0) - \oS° 2) 



O mixed 
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t,q,mixed 

t,q,mixed 
5 21 



-2Qf 1) + 3Qf 1) + Qf 0) + Qf 0) 



282 1 



g(001) _,_ ^' ^1 01)00 I 



91 
282" 



^(0000) 



Q (200) + Q (020) + A ( 4 ± v^3) Q(° 02) + £ (-5 ± V43) Q[ 



,(110) 



j(200) 



+ \ (l T 2^) QS 101) " ^ (17 ± 2v/43) u) - Q 
- ~ (4 ± V43") Qf 2) - g (-5 ± V43) Qf 0) - ~ (l T 2 Via) Q 



(011) 
2 



17 ± 2V43 Q 



.(101) 



19 =f 2V43 



Q (100) 



gCoio) 



+ 234 ( 33 T 16V ^) (^ f0 °° 0) " ^°° 00) ) ' 



(0000) 



^(0000) 
»2 



g,mixed 
2,0 



3 0} 



(0000) 



^(0000)\ ^ 



T-iq,mixed 
E 21 



79/9. 



(5.48) 



The corresponding eigenvalues (anomalous dimensions) are listed in Table |3J 

Closing this section we want to emphasize that all of the above results are written 
for generic quark flavors and can be applied to the study of baryon states with arbitrary 
quantum numbers. 



6 Nucleon Distribution Amplitudes 
6.1 The leading twist-3 distribution amplitude 



For completeness and for further use in Sec. |6.3.2| we write down the expansion for the 
leading twist-3 distribution amplitude $ 3 (x,/i) ( |3.2 ) 



$ 3 (X, /I) =XiX 2 X 3 C Nq 4>Nq{^) PNq{x) 
N,q 



(6.1) 
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where \i stands for the scale dependence. The polynomials Pjsr >q (x) are related to the 
eigenfunctions *&N,q of the Hamiltonian Hi/ 2 (4.49) and are given by 

PN, q (x u x 2 ,x 3 ) = (e^ XkZk \^! Nq {z)) , (6.2) 

where (...) is the 577(1, l)-invariant scalar product (4.34) corresponding to the spins 
ji — ]2 = J3 = 1- The normalization constants cjv, 9 are defined as 



-Nq 



|*Ar^(^)|| 2 /r(27V + 6) = Vxx^x^P, 



Nq | 



(6.3) 



Finally, the scale-dependent coefficients 4>N,q{p) are defined as reduced matrix elements 
of the multiplicatively renormalizable twist-3 operators 



'(/') =p Nq (d z )o tw - 3 (z^)y =0 , 



'Nq 



Nq- 3 (ri\N) = + \{pn)Nl (-ipn) N <f> Nq (») . 



Taking into account Eq. (6.3) one can project the coefficients 

<j>Nq{y) = J T>xP Nt q(x)^ 3 (x,fJ l ). 

The first few terms in the expansion (6.1 ) read 

i(2/3) 



as follows 



$3(^1,^2,^3) =120xix 2 x 3 
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(.(26/9) 

1,0 



P lfi (x) + U^J m P 1A (x 



63 

id 



,(38/9) 
J 2,0 



63 
~2 



+ ^4r P2 , l (x)+u^p 2 , 2 (x)+... 

5 " 



(6.4) 
(6.5) 

(6.6) 
(6.7) 



where 



2 
1 



(xi + x 3 - 2x s 



The superscript in 



7*2,0(3^) =3x\ — ?>x\x 2 + 2x\ — 6x1X3 — 3^2X3 + 3X3 , 
72,i(x) =(xi - x 3 )(xi + x 3 - 3x 2 ) , 
72,2(^) — x \ + x\ — \2x\X 3 + 9xiX 2 + 9x 2 x 3 — 6x3 . 
^Nq^ shows the corresponding anomalous dimension: 
a.(jJL) N En ^° 



(6.8) 



(6.9) 



These expressions agree with the ones existing in the literature, e.g. [HU [HI EH S3- In 
notations of Ref. [36] 

7 , fin/Sl ,j_ , „ , C2fi/9~l 21 

2 ( 



fc(2/3) 
J 



120! 2 o 6/9) 



,(10/3) 
1,1 



(10/3) 



(6.10) 
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6.2 Twist-4 distribution amplitudes: General formalism 

The construction of twist-4 distribution amplitudes is somewhat more cumbersome be- 



cause O(z) (4.6) involves five independent light-ray operators (for both chiral and mixed 
cases), and also because triangular mixing between quark and quark gluon operators 



makes the Hamiltonian non-hermitian. Solving Eq. (5.36) one obtains the expansion of 
O(z) in the form 



wherdU 



"Nq 



Nq 



J2 C Nq^Nq(z)Q Nq , 

N,q 



n-i 

U Nq 



(6.11) 



^Nq^Nq). 



(6.12) 



We remind that \Izjy are the eigenfunctions of the adjoint Hamiltonian (with respect to 



the scalar product (4.35)), cf. Eq. (5.37). It is convenient to split the sum in (6.11) in 



two: one contains the "three-quark operators", and the other one - three-quark-gluon 
operators: 



Nq[ 



"Nq 



+ B ^ 



Nq\ 



"Nq 



(6.13) 



N,q 



N>2,q 



Due to a block-triangular form of the Hamiltonian H, Eq. (4.5), only the first three 
("quark") components of the coefficient functions ^^(z), a = 1, 2, . . . , 5 corresponding 
to the "quark" operators are nonzero, ^ 



a=4,5 
N,q 



Nq\ 

0. On the other hand, all five components 



of the coefficient functions ^ 9 Nq are nonzero, in general. 



For the eigenfunctions of the adjoint operator W\ (5.37), the situation is the opposite: 



the "quark" eigenfunctions have all components nonzero, whereas "quark" components 



of the quark-gluon eigenfunctions vanish, \l/ 



g,t,a=l,2,3 
Nq 



0. Since the diagonal blocks of 
the Hamiltonian (4.5) are self-adjoint operators one can choose the eigenfunctions \1/ and 
as follows 



^Nq(z) 



Nq \ 





(6.14) 



Here N q (z) , ^ 9 Nq (z) are three-dimensional "vectors", whereas ^ 



two-dimensional. It follows from (6.13) and (6.12) that the only remaining "daggered" 



Nq\ 



Nq{z) are 



The coefficients Cf^q are added (cf. Eq. (4.9 1 ) to allow for arbitrary normalization for the functions 
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function v&jy (z) can be expressed as follows: 



(6.15) 



i.e. explicit construction of the adjoint Hamiltonian is in fact not necessary. 

Our task is to write the expansion of the nucleon distribution amplitudes in con- 
tributions of multiplicatively renormalizable operators. To this end it is convenient to 
introduce auxiliary amplitudes 



a \P) 



-{fi\)(-ipn) na N\_ m [ Vxe- i{pn) ^ XkZk O a (x) 



4 



a = 1,2,..., 5. (6.16) 



The r.h.s. of Eq. (6.16) involves the nucleon spinor N^(N^) for the chiral (mixed) 
operator, respectively. The factors (—ipn) na are introduced for later convenience: For 
the chiral operators we choose n = (0, 0, 0, 2, 2), and for the operators of mixed chirality 
n = (0,0,1,2,2). It will always be clear from the context which operator, chiral or 
mixed, is considered, so that we do not show it explicitly. 

The "standard" nucleon distribution amplitudes introduced in Sec. [3] are related to 
the amplitudes as follows: 



E±(x 1 ,x 2 ,x 3 ) = O w (x 1 ,x 2 ,x 3 ) , 

for the chiral case and 

^4(2:1, x 2 ,x 3 ) =O w (x 2 ,x 3 ,x 1 ) 
$4{xi,x 2 ,x 3 ) =0 {2) (x 1 ,x 3 ,x 2 ) 
D i (x 1 ,x 2 ,x 3 ) =C (3) (xi,x 3 ,x 2 ) 



^4^X1, x 2 , x 3 



x A ) = {4) ( 



x 1 ,x 2 ,x 3 ,x±) 



(6.17) 



, X 2 ,X 3 , X4) , X 3 , Xi, X4), 

$4(^1, x 2 ,x 3 , x A ) =C (5) (xi, x 3 , x 2 , x 4 ), 



(6.18) 



for mixed chirality. 

We define reduced matrix elements of the multiplicatively renormalizable local oper- 
ators as 



It follows from Eqs. (6.16) and (6.12) that 



>Nq 



VxJ2 P N q (x)n ab O b 



X 



(6.19) 



(6.20) 



ab 



Here 



^(s) = <e E *"**|tt&>- 



(6.21) 
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where (. . .) a is the SU(1, 1) invariant scalar product (4.34) corresponding to the confor- 



mal spins of the function We also define the functions 



(6.22) 



The scalar product for the functions ^Nq , ^Ng that corresponds to conformal operators 
(i.e. ^Ng{ z ) and W N (z) are shift-invariant polynomials) can be written as 

(*lr<,IW =c N jvx J2fc a (x)P^ a q (x)n ab P b Ngl (x) = c^P^ 

ab 

where cn = T(2N + 5) and the integration measure is defined as 



Ng>) , 



Hx{x) = x 2 x 3 



fl 2 (X) = XiX 3 , 



fiiix) = fi 5 (x) = ]^XiX 2 X 3 x\ 



fl 3 (x) 



X\X 2 the chiral case 
^XiX 2 xl the mixed case . 



Since the functions ^f^ q and W N form a bi-orthogonal system one obtains 

O a (x, fi) = fl a {x) ^ °Ng ^Nq(fj) PNq( X ) > 
Ng 



where 



C Ng - ( P Ng\ P Ng) ~ C N \/ C N . 



The scale dependence of the reduced matrix elements is given by 

<PNq{f*0) • 



We now specialize to the particular cases of interest 



(6.23) 



(6.24) 



(6.25) 



(6.26) 



(6.27) 



6.3 Chiral amplitudes S4, S| 

The expansion for the chiral quark, S 4 , and gluon, nucleon distribution amplitudes 
reads 



i(x,n) =x 2 x 3 



^2a Nq ^ Nq (fl)U Nq (x) + b Ng^Ng(^)^ 



Nq 



x) 



N,q 



N>2,q 



El(x,fl) =-XiX 2 X 3 xl Y b Nq^Nq(^)'^ 



N,q( X ) > 



(6.28) 



N>2,g 
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where U Nq (x) = P^= 1 (x), U. 9 Nq (x) = P$° = (x), U 9 Nq (x) = P^ =1 and the expansion co- 
efficients £jvfc, £ff k are defined as reduced matrix elements of the multiplicatively renor- 
malizable operators 

(0\O^ al+ (^)\N) = - l^X)m N Ni (-ipnffrM , 

(0|O^ + (/i)|iV) = - \(v\)m N Nl (-ipn) N e Nq (») ■ (6.29) 

Explicit expressions for the operators O c ^ al+ , O^J* rai+ for iV < 2 are given in Eq. (p.43h; 
the corresponding anomalous dimensions can be found in Table [2j 



Taking into account the symmetry properties (3.9), (3.10) one can write the normal- 
ization constants as 

a N \ =2 Jvx^(x) IL N>q {x) (2 + P 23 ) U N<q (x) , 

6^=2 fvxii4x)U^(xj(2 + P 23 ^U 9 Nq {x), (6.30) 

where P 23 is the permutation operator: P 23 $(xi, x 2 , x 3 ) = $(xi,x 3 ,x 2 ). Projecting out 
the contribution of a particular operator one finds 

\e Nq = J T>xU^(x)(2 + P 23 ) Ef (x) , 

\Zn 9 = J VxU Nq (x)(2 + P 23 ) H 4 (x) - § ' ( 6 ' 31 ) 



where 

m N 

" l qq 



, = 2 j Vx^{x)n Nq {x)(2 + P 2 ^n% ql {x). (6.32) 

Finally, taking into account the contributions with iV < 2 we obtain the distribution 
amplitudes 

28 ~ 

H 4 (x,/i) =4x 2 x 3 [eo,o(/i) n (x) + 9£ M (jjl) Hi (x) + 12&,o(p) + y£f, M n|(x) 

Sf (x, n) =^8! 0:1X2X3X4 £f (^)IIf (x) , (6.33) 
where 

n (x) =1, 

ili(x) =Xi + x 3 - -x 2 , 

n 2 (x) =x\ — 4xiX 2 + 2xg + 2xiX 3 — 4x 2 x 3 + x\ , 

—43 _ 13 

n|(x) = — x\ + 4xiX 2 — 2x2 — 47xix 3 + 4x 2 x 3 H X3 . (6.34) 
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and H| (x) = 1/2. In notation of Ref. [36] 

A2 =£0,0 , 

Ai/ 2 d =^o,o + ^i,o- (6.35) 
6.3.1 Mixed chirality amplitudes $ 4 , \I>4, 

The distribution amplitudes $4, _D 4 have collinear twist 4, but receive contributions 
both from geometric twist-3 and twist-4 operators. Moreover, among twist-4 operators 
there are descendants of twist-3 operators which matrix elements do not involve new 
nonperturbative parameters, cf. [46]. In the discussion of the operator renormalization 
in Sect. 5 such operators were excluded by imposing appropriate symmetry conditions 



on the solutions of the renormalization group equations, Eqs. (5.23), (5.32), and they 
are not listed in Sect. 5.4. They do contribute to the distribution amplitudes, however. 
The contributions to collinear twist-4 amplitudes that can be expressed in terms of the 
leading twist distribution are usually referred to as Wandzura-Wilczek contributions, 
with the remaining parts being "genuine" twist-4: 

$ 4 (x) =$f w {x) + $f -\x) , 
* 4 (a0 =K W (x) + ^T' 4 (x) , 

D 4 (x) =Df w (x) + Df-\x). (6.36) 
In what follows we consider these two types of contributions separately. 

6.3.2 Wandzura-Wilczek contributions $f w , #f w 

Let 

O tw ~ 3 (z, A) = O tw - 3 (z u z 2 , z 3 ; A) = O tw - 3 ( Zl , z 3 , z 2 , A) (6.37) 
be the leading-twist light-ray operator that enters the definition of the leading twist-3 



distribution amplitude: In comparison to (6.4) we have replaced z 2 <-> z 3 and added an 



argument A to remind that the "plus" projection is done with respect to this particular 
spinor. 

The short-distance expansion of O tw ~ 3 (z, A) can be written as 

O tw ~\z, \) = J2 C Nq U Nq {2) 0^7 3 (A) + 4 {n \ 3) S + *nM *[ P aa> ®X 3 W\) +■■■■ 

(6.38) 

Here the first term in parenthesis corresponds to the contribution of the conformal 



twist-3 operators (A) (6.4) which are annihilated by the step-down operator of 
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the SL{2) algebra, so that z[K M ^, 0^~ 3 (A)] = 0. The corresponding coefficient func- 
tions $N q (zi, Z21 z 3 ) are shift-invariant homogeneous polynomials of degree N and are 
related by a simple change of arguments Z2 «-> z 3 to ^N q (z) appearing in (6.2). The 
normalization coefficients are chosen as C Nq = \\§ Nq \\~ 2 = c Nq /T(2N + 6), cf. (|6.3|). 



The second term in (6.38) corresponds to contribution of the operators that in- 



clude one total derivative (in "plus" direction). Note that the corresponding coeffi- 
cient functions can be obtained by application of the three-particle "step-up" operator 
5^ = i + 5^=1 + ^=1 — J2k( z l^z k + 2^) (2.40) so they are not independent and 



S+ 



do not need to be calculated separately. The terms with two and more total derivatives 
have the similar structure. They are indicated by ellipses. 

Omitting contributions of genuine geometric twist-4 operators the expansion of the 



light-ray collinear twist-4 vector-operator 



t,mixedi 



z) reads to the same accuracy 



r(z, f i,\)] tw _ 3 =Y,c Nq y a Nq (z)o Nq ( f i,\)+ c N,n 

1 



Z) U Nq {}X, Aj + 2^ Un 9 ^N +ltq (z)0 N+1>q (fM, A) 

N,q N>l,q 



Nq 



N,q 



2(2iV + 5)' 



S + V a N (z) z[P xx , Nq (v, A)] + . . . , (6.39) 



where 



v Nq 



»N+l,q 



CM) 



1 



(M)o^- 3 (A) 



N + 2 



1 



N + 2 



XX > ^Nq 



CM)] (6-40) 



are the light-cone conformal operators ifK^, Qjv (/-f, A)] = ifK^, Qjv+i, (A, fi)] = 0. (It 



is tacitly assumed that the generator S + in the second line in Eq. (6.39) involves the 
conformal spins of the functions it acts on.) The operator Ojv ? (yU, A) has geometric 
twist-3 whereas Ojv+i, g (/i, A) contains both twist-3 and twist-4 contributions. 

The task is to find the corresponding (five-component) coefficient functions ^Nq, ^Nq- 
This can be done by observing that O a (z, fi, A) only depends linearly on \i and essentially 
reduces to O tw ~ 3 (z, A) after a formal substitution // — > A. Taking into account the 
identities 

\d^O {a=1 ' 2) {z,n,\) —O tw ~ 3 (z, A) , Xd^O {a=3 \z,fi\) = -d Z2 O tw ' 3 (z, A) (6.41) 



and comparing the two representations in ( ]6.38[ ) and ( J6.39[ ) one easily gets for ^^q 

(6.42) 



--®Nq{z) 



o=3/ 
Nq \ 



'3^Nq\ 



,a=A 1 



f a=5, 
Nq I 



0. 



For the functions ^% q one obtains, after some algebra 



^ N+l,q (ZJ 



S + -2(N + 3)z 1{2) $ Nq (z) 



2(2N + 5) - [S+ - 2(N + 3)z 3 ] d 3 



Nq[ 



(6.43) 
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where S + = S^ =1 + S^ 2=1 + 5^ =1 . One can easily check that the functions defined 



in Eq. (6.43) are shift-invariant polynomials, as they should be. They are normalized as 
\\V Nq \\ 2 = {N + 2)\\$ Nq \\ 2 , \\$n+i*\\~ 



2(2A^ + 5)(iV + 3)'||$ 



Nq\ 



(6.44) 



Going over to the "P-representation" as defined in Eq. (6.22) we obtain 

PSlW =dxx 1 P^{x) , P$(x) =d 2 x 2 P^7 3 (x) , P$(z) = - 2d 3 P^ q -\x) (6.45) 
and 



^N+l,q( X ) 



2N + 5 — Xd Xl J xi P, 



ytw— 3 1 
Nq 1 



X) 



■■{2N + 5- Xd x ^jx 2 P^~\x) 



2(2N + 5-Xd 3 )p t Z- 3 



[x 



(6.46) 



J 3J i Nq 

where X = x\ + x 2 + x 3 . 

Taking appropriate matrix elements we end up with the Wandzura-Wilczek contri- 



butions to the auxiliary distributions (6.16) 

CNq 4>Nq 



WW 



/Ji a yX 



E 

Nq 



P\ r „ ix) 



(27V + 5) \N + 2 Nq 



P 



(a) 



N + 3 N+1 ' q 



x 



(6.47) 



where CN q ,(pNq are the twist-3 expansion coefficients, (6.1). One can easily verify the 
following relation (cf. Eq. (3.19)): 



[0^{x)\ 



WW 



X\ 



pW(x)] ww + x 2 [0^(x)\ 



WW 



Finally, taking into account Eqs. (|6.18|) we obtain 

CNq 4>Nq 



WW, 



X) 



-E 



*f w (x) 



N,q 



E 



{N + 2){N + 3) 

CNq 4>Nq 



(N + 2)(N + 3) 

ytw— 3 
Nq 

one of the main results of this paper. 



N + 2- 



N + 2- 



d 

dx 3 

d 

dx 2 



XlX 2 X 3 Pjfg "{Xi,X 2 ,X 3 ) 



(6.48) 



(6.49) 



xix 2 x 3 P^ q 3 (x 2 ,x 1 ,x 3 ) , (6.50) 



N,q 

where the polynomials -P^~ 3 (x) are defined in Eq. (6.1 ), (6.8). These expressions present 



6.3.3 Genuine twist-4 contributions " 4 , $f" 4 , tff, $1 



The expansion of the twist-4 auxiliary amplitudes O a (x) (6.16) reads 
O a (x) =ti a (x) 



^A N qVNqV Nq (x)+ BNqVNqKqi 
. Nq N>2,q 



X) 



a =1,2,3, 



O a+6 {x) =-XiX 2 X 3 X% 



B NqVNq^Nq( X ) 



a = 1,2. 



(6.51) 



N>2,q 
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The coefficients An iQ and q are given by 



A Nq 



1 

Xi 



V l Nq {x)\ 2 + -\V 2 Nn {x)\ 2 + ^\P% n (x 



x 2 



Nq 



Nq\ 



b n\ = I Vxfi 4 (x) ^V 9 / q {x) u ab V 9 / q (x) , 

. Note a useful identity 



(6.52) 



where in the last line u 



2 1 
1 2 



/•X3 

/ drr (9iP^ g (xi,x 2 ,r) + d 2 V^ q (x 1 , x 2 ,t)) (6.53) 
./o 



which is a consequence of Eq. (5.32). 

The reduced matrix elements r] Ng (fi) and f] 9 Ng (fi) are defined similar to (6.29), replac- 
ing the spin-down spinor A/4 by the spin-up one, N* , cf. (6.16), (6.19). The expressions 
for the operators ®™ q ed , Q,9^ xed f or jy < 2 and the corresponding anomalous dimen- 



sions are given in Eq. (5.48) and Table |3j respectively. They can be isolated by the 
projection 

3 

0=1 ^ g' 



(6.54) 



where 



M ^, = B Nq > I Vxxix 2 x 3 



qq 



(^Nq^nl^) + ^V% q {x)Vf ql {x) + ^V 3 Nq (x)Vf q ,(x)) . 

\X\ X 2 4 J 
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Taking into account the contributions of local operators with N < 2 we obtain 

r 20 

^ 4 w ~ 4 (x,/i) = 12x1X3 rjofiifj) + 4r} lfl (ii)Vi,o{x 2 ,x 3 ,x 1 ) + — 771,1 (/x) Vi,i(x 2 , x 3 , xi) 



+ 



+ 



5 /ll 5 

+ 



2 V 2 x/43 



45 

V2,0 W V 2fi (x 2 , X 3 ,Xi) + —r]2,l{lj) V 2 ,l{x 2 , X3,Xi) 



5 /ll 



2 V 2 v/43 
70 



^2, 2 (/i) ^2,2(2:2, &3, Zi) + ^?72,o(^) ^2,0(^2, ^3, &l) 



4 (x,/j) = -12xiX 2 ?7b,o(A*) +4771,0(^)^1,0(33,^1,^2) - 71/1,1(11)^1(13,11,12) 



20 



where 



+ 



+ 



5 /ll 



2 V 2 
5 /ll 



+ 



45 



. V2,o(fi) V 2 , (x 3 ,x 1: x 2 ) - —rj 2jl (fi)V 2 ,i(x 3: x 1: x 2 ) 
43 / ^ 



2 V 2 

70 
47 
1 



^ V2M ^2,2(2:3, X U X 2 ) + ^^f.oM ^O^, X ± , X 2 ) 



vL(i^)'PL(x3,x 1 ,x 2 ) 



^l(x,fx)— -8!xix 2 x 3 x 4 



$l(x,fA) = - -8!xix 2 x 3 x 4 



V2M + oV2,l(^) 



vl,M - 3^2,1 (^) 



(6.55) 



Pi,o(a:) =x x + x 2 - -x 3 , 
Vi,i(x) =Xi -x 2 + ^x 3 , 

2 

V 2 ,i(x) =x\ -x\- 2xix 3 + 3x 2 x 3 - -x\ , 



^2,0(3) 
^2,2 (x) 



~-x\ + - (-5 ± v 7 ^) x x x 2 +x^ + -(lT 2^43) X1X3 
-i(l7±2v / 43)x 2 x 3 + ^(4±v / 43) x 2 3 , 



~ 11 17 

V 2 (x) =64x1 - 55xiX 2 + — x\ - 73x1X3 + llx 2 x 3 + — x\ , 

V 20 (x) =16x^ x\ — 32xix 3 + x 2 x 3 + 5xg . 

3 



(6.56) 
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The expression for D±{x,i£) (3.15) can be obtained using the integral representation in 



Eq. (3.19) 



The reduced matrix elements r]^ q for N = 0, 1 are related to the parameters intro- 
duced in Ref. 1361 as 



Voo 



13 11 

Ai f x = - -0oo - tt^oo - rVw + q^n 
o 10 5 3 



Ax/r 



6 



^oo 



10 



Voo 



(6.57) 



The term in O o is the Wandzura-Wilczek contribution that can be traced to the twist-4 
operator containg a transverse derivative of a local twist-3 three-quark operator (and 
properly symmetrized). 



7 Conclusions 

The motivation for our study has been to work out efficient techniques for a calculation 
of anomalous dimensions of generic higher twist operators in QCD. Apart from the ap- 
plications to QCD phenomenology, this project was fuelled by the recent progress in the 
understanding of the spectrum of the dilatation operator in the maximally supersymmet- 
ric N = 4 Yang-Mills theory [26] and, in particular, the work [27] where it was argued 
that the diagonal part of one-loop QCD RG equations (for arbitrary twist) can be writ- 
ten in a Hamiltonian form in terms of quadratic Casimir operators of the full conformal 
group 5*0(4,2). In simple words, the symmetry under the conformal transformations 
in the directions orthogonal to the light-cone plane implies existence of relations be- 
tween the renormalization group equations of different (geometric) twist. This, in turn, 
suggests that the techniques developed for the description of quasipartonic operators in 
QCD [23] can be generalized to include non-quasipartonic operators as well. Our goal 
was to develop a consistent computational framework based on these ideas. 

The first step was to construct the appropriate operator basis with "good" transfor- 
mation properties. We found that the complete basis of one-particle conformal operators 
for chiral quark and self-dual gluon fields in QCD contains seven light-ray fields, and sim- 
ilar in the anti-chiral sector. An interesting feature of this basis is that it includes some, 
but not all, transverse derivatives: If the transverse plane is parameterized in terms of a 
single complex variable as it is usually done in the studies of high-energy scattering, then 
the basis fields only include holomorphic derivatives acting on holomorphic components 
of the fields, and vice versa. 

Although much of the formalism appears to be general, in this paper we concentrate 
on the simplest example of non-quasipartonic twist-four baryon operators that contain 
two "plus" and one "minus" quark field, schematically q + q_q + , and their mixing with 
(quasipartonic) four-particle operators involving a gluon field, of the type q + q + q + F + ±. 
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For this setup we calculate all one-loop evolution kernels and check that they are SL(2) 
invariant, as expected. The evolution equation for three-quark operators of the same 
chirality turns out to be completely integrable. The spectrum of anomalous dimensions 
coincides in this case with the energy spectrum of the twist-4 subsector of the SU(2, 2) 
Heisenberg spin chain, confirming the prediction of [27] . We find the explicit expression 
of the corresponding conserved charge and calculate its spectrum. A simple analytic 
expression is found for the lowest anomalous dimension for chiral quark twist-4 operators 
with odd number N = 2k + 1 of covariant derivatives. For other cases the spectra are 
studied numerically, see Figs. [2}|6] and Tables 2|3 It turns out that differences between 



twist-4 and twist-3 operators and also between twist-4 operators of different chirality 
mostly affect a few lowest eigenstates (for a given N); the upper part of the spectrum 
of anomalous dimensions is universal. The spectrum of twist-4 quark operators overlaps 
strongly with that of quark-gluon operators, apart from a few lowest states. 

Finally, these results are applied to give a general classification and calculate the 
scale dependence of subleading twist-4 nucleon distribution amplitudes that are relevant 
for hard exclusive reactions involving a helicity flip. In particular we introduce new four- 
particle distribution amplitudes involving a gluon field, and derive explicit expressions for 
the expansion of all distribution amplitudes in contributions of multiplicatively renormal- 
izable operators taking into account first three orders in the conformal spin expansion. 
As a byproduct of our analysis, we give an all-order expression (in conformal spin) for 
the contributions of geometric twist-3 operators to the (light-cone) twist-4 nucleon dis- 
tribution amplitudes, which are usually referred to as Wandzura-Wilczek contributions. 
The applications of these results to phenomenology of hard exclusive reactions will be 
considered elsewhere. 

The techniques suggested in this paper can have a rather broad field of applica- 
tions, in particular to the calculation of twist-4 corrections to the structure functions of 
deep-inelastic lepton-hadron scattering. We plan to consider this problem in a separate 
publication. 
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